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Prefac e to  . Chapters  3 and  _4 


Chapters  3 and  k complete  Part  I of  this  report. 

YJe  would  like  to  call  attention  to  the  fact  that 
the  informal  nature,  of  these  unpublished  Technical  Reports 
invites  an  amount  of  speculation  which  would  be  perhaps  in- 
tolerable in  the  published  literature,  and  a certain  loose- 
ness and  crudeness  in  derivations  and  formulations  of  prob- 
lems which  a regular  published  work  would  not  exhibit.  It 
seems  desirable  to  the  authors  to  make  a statement  at  this 
point  concerning  the  scope  of  what  we  know  and  what  we  do 
not  know  in  the  subject  matter  of  these  two  chapters. 

The  reader  will  quickly  see  that  the  subject 
matter  of  Chapter  3 is  confined  to  the  hydraulics  of  sharply 
stratified  media,  whereas  real  estuaries  are  always  more  or 
less  diffusely  stratified.  What  is  more,  no  discussion  is 

V-  **  ' * v.-se-  - -•*  — 

made  of  the  order  of  magnitude -of  the  friction  terms.  . Ip. 

no’T'^r  «a  i ■ n rrl  a 1 a Tr  a r»  *P1  rma  ( an  n h pa  -fn  r»  ^ vr  o r»  a j onrri  ngoT»g  a "1  — 

ready  have  crude  approximations  of  the  friction  terms  (Chezy 
and  Hanning  formulas),  but  we  do  not  have  even  these  approxi- 
mations for  two  layer  flow.  For  this  reason  the  differential 
equations  of  gradually  varied  flow  of  two  layers  are  for  the 
most  part  left  unintegrated  and  all  that  is  demonstrated  is 
the  qualitative  aspects  of  the  flow. 

In  the  case  of  entrainment  of  water  from  one  layer 
into  another  we  can  only  perform  integrations  of  the  equations 
when  the  amount  of  entrainment  is  known,  whereas  in  real  estu- 
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a-ries  we  do  not  have  a priori  Knowledge  nf  this  amount. 

The  reader  will  see,  therefore,  that  the  subject  flatter 
of  Chapter  3 is  really  very  incomplete,  leaving  undeter- 
mined all  the  constants  which  depend  upon  turbulent  mix- 
ing, upon  the  frictional  stresses  on  the  bottom, -and  the 
free  surface  and  the  walls,  and  upon  the  amount  of  entrain- 
ment. 

The  contents  of  Chapter  4 are' somewhat  iiffef- 
ent.  First  Of  all,  they . contain  summaries  of  several 
already,  published  papers  on  the  mixing  in  estuaries.  Most- 
of  these  papers  have  proceeded"  on  the'  basis  of  hypotheses 
about  the  nature  of  the  mixing  process.  The  applicability 
of  these  hypotheses  appears  to  be  restricted  to  only  cer- 

estuaries  , and  it  must  be  admitted  that  more  worK  has^  • 
been  do  he’- 'that  involves  guessing  what . the  mixing. processes  . 
in  an  estuary  might  be,  than  has  been  done  in  trying  to 
•find  out  what  the  mixing  processes  in  an  estuary  actually 

As  incomplete  as  the  subject  matter  of  Chapter  4. 
is,  it  is  hoped  that  it  will  suggest  which  of  the  possible 
mixing  processes  in  estuaries  may  be  important  in  any  par- 
ticular one  which  is  -the  subject  of  study,  and  that  it  will 
also  suggest  the  typo  of  observations  which  will  be  most 
desirable  in  studying  a particular  estuary.  For  example: 
in  an  unstratified  estuary  it  seems  that-  a more  or  less 
uniform  spacing  of  stations  u.p  and  down  the  estuary  is  de- 
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sirabie;  but  in  an  estuary  vvliich  appears  1 d'  be  "subject  to 
'ttL&r Constraint  of  overmlxing  "(Section  4.51)  the  location  of 
stations  shouid.be  largely  confined  to  control  sections. 
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CHAPTER  3. 

Gradually  varied  flow. 

3.1  The  general  equation  of  varied  flow. 

We  envisage  a channel  whose  axis  is  oriented  along 
the  2L  axi9>  whose  width  is  ).  The  axis  i3  directed 

vertically  upward.  The  ^ - axis  is  directed  horizontally 
across  the  channel. 

Two  layers  of  vertically  homogeneous  liquids  flow 
in  this  channel,  the  lighter  on  top  of  the  heavier.  The 
upper  liquid,  whose  density  i3  j?  ( y,  ),  has  a free  surface  at 

2-  = S,  (*). 

The  interface  between  thi3  upper  fluid  and  the  low- 
er liquid  layer,  whose  density  i3  ( *X  ) , is  at  * - )• 

The  bottom  of  the  channel  is  at  y£  ~ (?*)• 


Z 


* 
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For  convenience  auxiliary  quantities  are  intro- 
duced : 

the  layer  depths: 

Dl  = 5,  - 4u  D2  = '■** 

the  discharges  per  unit  width: 

°-l  = Dlul  ^2  = D2u2 

and  the  total  discharges: 

Ol  = bQi  0.2  = ^2 


The  steady  state  equation  of  motion  in  the  ^—dir- 


ection is 


-t 


2. 

a-t  / a* 


The  viscous  term  is  supposed  to  be  primarily  ver- 
tical (as  might,  be  the  case  in  a wide  channel  with  thin 
layers)  and  it  is  supposed  that  A4"  the  dynamic  eddy  viscos- 
ity,  may  be  a function  of  £ . Therefore  the  term 
is  used  instead  of  the  usual  yu.  VV/U- 

The  additional  term  — (<L)  is  intended  to  in- 

clude retarding  forces  such  as  might  be  due  to  a fine  screen, 
or  pilings,  or  long  grass. 

The  steady  state  equation  of  continuity  is 


(2) 
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Now  we  know  that  the  following  identities  hold: 


>U. 

- (W 

(3) 

p/U“  

- ^ . £( O'lL'ir) 

^ " 

- ~ ^ ^ w 
N . 

(4) 

' 

QUJ-  - 

J 

- f p 

**  r 

{ 5) 

By  adding  these  three  the  last  terms  vanish  on  account  of 
the  continuity  relation,  and  we  may  substitute  the  remain- 
ing terms  of  the  right-hand  member  into  the  left-hand  member 
of  the  equation  of  motion,  obtaining  the  following  form: 


We  now  integrate  this  equation  vertically  over  each  layer. 
First,  over  the  top  layer: 


(7) 
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The  limits  and  X are  functions  of  . 

Differentiation  of  the  following  integral  yields  the  follow- 
ing: 

2.1, 


f. 


(C) 


The  vertically  integration  equation  of  motion  becomes: 


Also,  we  may  integrate  over  the  width  of  the  channel  hor- 
izontally along  T-  and  dividing  the  result  A obtain 

Cr  * 


the  following  expression: 


3 


r 


& 

St 


> I 


(r^l 


i 


fv 


>?,  110) 
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At  the  solid  side  walls 


/ ^4- 

2- 


Therefore  the  above  equation  of  motion  becomes 


br 


n 


(r+r  ) 


. -L^r  ?,  P,  ^ 1- 


I 


I 


*£&.*■>?)  I 

dr  J 


in 


_f.  p/W.  t-cT 


a* 


ties  I 

- (r^+r)jg 


r** 


K 


(11) 


We  now  study  each  of  the  two  layers  separately. 
Suppose  that  within  the  top  layer  is  uniform  with  ^ , 
that  is,  has  a value  v;e  may  designate  as  . The  hydro- 
static pressure  is 

?•  = W(r' 

Therefore 


(12) 
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If  no  fluid  comes  through  the  top  boundary 


W 


**  *, 


(13) 


or  more  simply 


© 

s. 


f 


w 


<SK 


J= 


-atv 


s\ 


(14) 


At  the  interface  3^-  ~ c *)  there  may  be  some  flow 
of  deep  water  up  into  the  upper  layer,  so  that  if  CUZ, 
is  the  upward  velocity  of  water  penetrating  the  inter- 
face we  can  write  our  expression  as 


(15) 


The  integrated  equation  of  motion  of  the  top 
layer  then  becomes: 


^ - r,  -tl  - / 


(16) 


where 


r.-  -/“■** 


^ 2--/, 
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In 

J 


ir% 


the  bottom  layer  * 

A,  * * /-  O 

O 


*li 

61 


(f~'+P)j, 


Sc 


f, 

4 


This  expression  becomes 


<0  V/" a A +#x  r/ 


+ f- ^ 

, j_  d£z  ?•>«'*  ^ ^ Tv “ “ 

f *fr-  **  J 


(17) 
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3. 2 A single  fluid:  effoot  of  friction 

Let  us  first  consider  a single  layer  of  water  of 
uniform  density  ^ , and  consider  the  surface  as  the 

solid  bottom.  Let  us  suppose  that  there  is  no  frictional 
stress  at  the  free  surface  ( T,  * 0 ) and  that  there  is  no 
addition  of  water  to  the  flow  (£»■>£*.  *o).  The  equation  (1) 
then  becomes 


£ (a  t . *.*• 


,2>, 


(i) 


= — T- 


Furthermore,  let  us  first  consider  a simple  chan- 
nel of  constant  width  ( ).  By  continuity 

. The  equation  may  be  written 


or 


#(-£***J 


(2) 


The  condition  of  flow  in  which 


/ olf. 


is 


called  "uniform  flow"  by  hydraulic  engineers  (Bakhmeteff, 
1932).  According  to  the  Cnezy  formula  for  very  wide  channels 

'Z'C;  * C.r^~Sx,  ^ ~ 

or 

- D,  CJTL 


l 
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where  C is  the  Cnezy  resistance  factor  determined  by  empir- 
ical formulas.  We  designate  a quantity  IK  -2,C  as  the 
conveyance 


= 6' 


n) 


This  corresponds  to  the  case  in  our  formula  where 


= 2- 


These  two  formulas  can  be  reconciled  if  we  write  "Z\  * 
thus 

K = f ~&t  s*- 


or 


K'U., 


=.  5*- 


T ot 


or  no  on 


1 Tr  o 


ot-U,  


JK  * 


D, 


-V/, 


(4) 


We  thus  see  that  the  depth  of  the  layer  increases 


or  d6cr6°£63  depend inr  upon  ths  si.^n  q£  ^otii  numomt/OT*  o.nd[ 
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denominator.  The  case  where  the  denominator  vanishes  is  a 
critical  one;  it  occurs  at  a critical  velocity  of 


the  velocity  of  propagation  of  a long  gravity  wave.  If  the 
stream  is  flowing  with  a velocity  less  than  the  critical  one, 
(i.e.  , tt  < ) the  flow  is  subcritical  and  the  denominator 

is  positive.  For  supercritical  flow  the  denominator  is  nega- 
tive. It  is  interesting  to  note  that  for  subcritical  flow  in 
a channel  with  level  bottom  { s ° ) the  effect  of  friction 
on  the  Dottom  is  to  decrease  in  the  downstream  direction 

and  hence  to  accelerate  the  flow.  The  gradually  varied  flow 
in  most  natural  v/a tercourses  appears  to  be  subcritical,  and 
we  will  confine  our  attention  to  such  regimes. 
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! 


3. 3 Single  layer:  widening  the  channel 

Let  us  now  consider  the  effect  of  widening  of  the 
channel  alone  on  the  flow.  The  equation  (3.1-16)  becomes 
simply 


di-  ( ?>, 


I 46-  t>'~ 

+ -3T 


= o 


(1) 


Now  'ey  continuity 


=.  6 


(2) 


Differentiation  of  the  equation  of  motion  yields 


•LV, 


i 

4r  df.  0), 


(3) 


The  result,  in  subcritical  flow,  of  a broadening  in  the 
channel,  is  to  thicken  the  single  layer. 


Single  layer:  entrainment 


age 


3. 


If  the  only  effect,  permittee  is  that  mass  be 
added  uo  the  single  layer,  (i.e.  O ) we  obtain 

the  following  equation: 


This  corresponds  to  a case  where  water  of  density  13 
injected  into  the  single  layer  with  a velocity  of  and 
per  unit  width. 

By  continuity 


U 


i a 


(2) 


(3) 


This  is  approximately 


fki 

**- 


r 


1 *~§'U4. 


U) 


Now  going  back  to  the  eouation  of  motion 


21>  - AjI—  ( f.  °A.  -+-3),  i£i) 

2>,f,  e^/L  Cb.f,)'-  r ‘2”' 


(5) 
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First  consider  the  case  where  , then 

-^-£r+&lVl  Hr  -~fl  ^ ^ 

T>,  J),  f,  ** 


(6) 


The  equation  Tor  the  change  of  with  % is: 


etD,  _ p,  ^ ^ 


The  result  is  that  in  subcritical  flow,  as  long 
as  the  injected  water  is  injected  with  a horizontal  velocity, 
nt-i.  , more  than  the  stre^-n^t*.  , the  layer  will  thicken. 

The  next  step  is  to  compare  these  effects  in  some 
special  cases  where  a bottom  layer  is  present. 


I 


* 


I 


t 


\ 
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3. 5 Two  layers:  deep  bottom  layer 

Suppose  the  ueep  layer  does  not  move  at  all,  and 
that  it  does  not  .nix  with  the  upper  layer.  Ir.  this  case, 
equation  (3.1-17)  takes  on  a particularly  simple  form: 


+■  (f.V,  4 /h.  Vl) 

Thus 


D,* 

iy 

tk 


(l) 


a-  o 

(2) 


i 

I 


or 


cL 

aOj/L 


(fV) 


(f^  &) 


o 


(3) 


This  equation  essentially  states  that  the  horizontal 
pressure  gradient  in  the  bottom  layer  vanishes. 

We  may  now  combine  this  expression  with  the 
equation  (3.1-16)  to  obtain 


-t 


-L 


- - Vi 


(4) 


1 


or  introducing 
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i 


A?tf 


oLV\  _ 

cAsjL 


f/L  _±_c*6- 

- +■  -r—L  Zl7 

oxp,1*  ~ ^ 


(5) 


D - til 


The  denominator  defines  a new  critical  velocity 


*14, i VLP;  which  is  considerably  less  than  that  for 

one  layer  flow,  su-c  = vpT  . Velocities  of  flow  whose 

rr  tc, x 

interfacial  Froude  number  r i s — ■■ ■ is  less  than  unity 

£ 

we  may  call  subcritical.  It  is  seen  that  the  same  general 
effects  that  occur  in  the  subcritical,  one-layer  case,  Fj  < I 
occur  in  subcritical  (referred  to  ) two -layered  flow, 


except  that  the  behavior  of  the  top  layer  is,  so  to  speak, 
inverted  vertically. 


Section 


3.  ^'i  { rage  i 


3, 51  Effect  of  entrainment 

Suppose  that  the  lower  lc-yer  is  very  deep,  but  that 
mixing  occurs  across  the  interface.  The  model  therefore  is 
of  this  kind: 

(i)  The  interface  g-  = permits  upward 

movement  of  deep  water  into  the  top  layer  where  it  is  mixed, 
but  no  mixing  downward  is  permitted.  Thus  > Pi  13  indepen- 
dent of  ft  , but  Ot  'j  as  oO  . The  mix- 
ing may  be  due  to  winds,  tidal  currents,  or  the  shear  devel- 
oped at  the  interface,  but  in  this  model  it  is  regarded  as 

fixed  independently  of  the  mean  flow. 

(ii)  The  depth  of  the  top  layer,  d,  is  very 

much  less  than  that  of  the  bottom  layer. 

(iii)  Y/ithin  the  top  layer  mixing  is  so  strong 

vertically  that  the  density  the  velocity 

are  independent  of  . 

[nv>  A IfAl  AA<  ffT  A 'P  V •?  n r»  A ^ A A A 1HO  f O ^ n t A f Vl  I I r.  AQ  T» 

•A  AAV/  » WX  WW  A WJ  W A UtA  /W  A A A W A UOV  y KVA  WA  -A.  A - W V/  w AA  V VA  v/  A 

layer  is  UK**  . The  mass  flux  per  unit  width  of  the  top 
layer,  ID,j5y/U/  , increases  by  the  addition  of  deep  water 
at  the  rate  so  that  by  mass  continuity  the  following 

relation  must  obtain: 


0 

(fit 


(&?>«■•)  = 


^ ft 


(I) 
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If  the  salinity  of  the  two  layers  are  denoted  by  .5 
and  Sj.  * respectively,  the  salt  transfer  per  unit  width 
in  the  top  layer , J),  u»  5,  increases  by  addition  of  salt  water 
from  below,  and  by  conservation  of  salt  the  following  rela- 
tion is  obtained: 


eW 


'UTZx 


(2) 


If  the  density  is  simply  a linear  function  of  salinity, 
f 1 equations  (1)  and  (2)  combine  to  yield 

the  following  form  of  relation: 


57  hW  " 


(3) 


Because  of  its  great  depth,  the  velocities  in 
the  deep  layer  are  small,  and  hence  the  horizontal  pressure 
gradient  at  any  depth  within  the  deep  layer  must  vanish. 
This  is  expressed  by  an  earlier  equation: 


d=-  (f,  V,  = 6 


U) 


* ° t i O’"*  ^ ^ 1 ( ^“>-0  r-o 


Equation  ( i J simplifies  to  the  following  form: 


^ ^ L^r-]  - >f'u 


cL^. 

15) 


W 

T3  tr  olnm^nof^rv*  r\  ■f  S V>o  Vuto  on  onMr’t’^ona  ( I ^ u ^ C \ ♦-  r* 

— v‘  - ^ - 1 \ * - **•■*  v x * 


following  equation  is  obtained 


J- 

tty 


2\?< 


Y 


= o 


(6) 


where 


■i  = 


The  meaning  of  equation  (6)  can  be  explained 
easily  if  we  introduce 

A "H  A O t 

st\  - r* 


(V) 


a function  denoting  the  transport  of  the  upper  layer. 


Clearly 


•v  » 


(3) 


where  the  subscript  q indicates  values  of  Y and  ^ at  <XsO 


v O 
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The  equation  (7)  then  becomes 

jLFJl  + illl  „ 

^ p J 

or  if  K — ^ ~Y'A  ^ 


(9) 


(10) 


Changing  variables 


c tLD  _ 


'W.l 


2^-« 

/it* -2* 


(id 


or  if 

- '«•*/ K 


(12) 


g4  D, 


2-6- -I 

F^- 


I 1 7) 
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The  numerical  values  tabulated  in  Table  3.51.1  in- 
dicate that  there  is  a range  of  0,  f < Q , for 

which  the  depth  of  tue  top  layer  decreases  with  increase  of 
transport.  In  the  case  of  Alberni  Inlet  this  is  apparently 
what  happens.  hloreover,  the  value  of  increases  rapidly 
after  a certain  value  of  the  transport,  so  that  evidently 
the  two  layer  model  breaks  down  where  ^ ?.0 

The  veloci  ty,  C.  , of  an  internal  wave  at  the  in- 
terface is  given  by 


9 


vD, 


From  equation  (12)  one  sees  that  the  velocity  of 
the  upper  layer,  nJ-t , is  equal  to  C at  the  value  Z . At 
points  further  upstream  C <• 

Tully  has  communicated  to  the  writer  the  fact 
that  the  internal  wave  motions  of  the  open  sea  do  not  seem 
to  p,0noti*s.ts  n tiis  two  nort ion  q ^ a Hoop  ysred 

estuary.  Apparently  this  point  of  critical  velocity  acts 
as  a block  to  such  deep  ocean  waves  and  prevents  then  from 
progressing  into  the  estuary. 

This  leads  to  an  interesting  result.  The  mouth 
of  the  estuary  will  act  as  a control  so  that^-L  at  the 
mouth.  The  depth  of  the  upper  layer  is  a maximum  at  a 
this  corresponds  to  the  place  where  the  velocity  of  flow  is 
one  half  the  critical  velocity. 


TABL2  3.51.1 


b 

0 


u 


2b-l 

b-2 

0.  5 


. 2 


0.3 


• u 
. 5 
.6 
.8 
1.0 
1.2 
1.4 
1.6 
1.8 
2.0 
2.2 
3.0 
10.0 
100.  0 


0.  125 
0.  0 
-0.25 
-0.  50 
-1.  00 
-1.75 
-3.00 
-5.50 

-13.00 

- 

+ <00 

+13.0 

+ 5.0 
2.4 
2. 01 
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3. 52  Experimental  studies  of  entrainment  In  a two  layer  system 

In  the  experimental  flume  it  is  possible  to  add  either 
fresh  or  salt  water  to  the  top  layer  by  means  cf  a shower  bath 
falling  onto  the  free  surface.  This  is  an  artifice  for  increas- 
ing the  discharge  of  the  top  layer  along  the  axis  of  flow.  It 
differs  from  the  natural  entrainment  process  ir.  estuaries  in 
two  ways:  (1)  In  estuaries  only  salt  water  is  available  for  en- 
trainment, but  experimentally  we  can  add  either  fresh  or  salt 
water;  (2)  In  estuaries  a counter-flow  occurs  ir.  the  deep  water 
due  to  the  entrainment,  whereas  in  the  flume  the  water  is  added 
externally,  so  that  it  does  not  cause  a reverse  flow  in  the  deep 
water,  but  such  a flow  can  be  induced  by  an  auxiliary  pump. 

Equation  3~sr.u  may  be  written  in  two  forms  depending 
on  whether  the  entrained  water  is  fresh  or  salt, 
fresh  water  entrained 


Salt  water  entrained 


(2) 
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These  two  quantities,  ? (K)  and  Jare  given 

in  Table  3.52.1.  This  table  is  perhaps  somewhat  easier  to  use 
in  computing. 

The  results  of  experimental  runs  in  the  flume  are 
given  in  Table  3.52.2. 


/ 
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R 

0 

0.1 
0.2 
0.25 
0.3 
0.4 
0.  5 
0.6 
0.7 
0.8 
0.9 


Table  3.52.1 

?(F) 

a(F) 

0.  00 

-0.22 

+0.33 

-0.50 

+0,13 

-0. 66 

0.00 

-0.68 

-0.15 

-1.  36 

-0. 50 

-2.0 

-1.00 

-3.0 

-1.74 

-4.6 

-2.95 

-8.0 

-5.50 

-18.0 

-13.00 



— oO 

1.0 
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Table  3 

.52.2 

Run 

« 

I 

1 

2 

3 

4 

5 

River  discharge  per 
unit  width  cm^sec"! 

15.8 

31.4 

31.4 

11.4 

i 

11.4 

Shower  discharge  per 
unit  width  cm^sec"! 

7.9 

12.8 

12.2 

11.4 

11.4 

i 

Shower  Fresh  or  Salt 

S 

F 

S 

F 

^ I 

9.0 

10.0 

10.  0 

14.2 

1 

13.8 

V 

8.0 

8.0 

S.5 

13.8 

16.2 

B1 

720 

1000 

1000 

2650 

2620 

B.'3 

510 

510 

610 

2640 

4250 

/■ 

1.009 

1.005 

1.  004 

1.009 

1.009 

/ 

1.018 

1.011 

1.014 

1.009 

1.016 

•/ 

V ' 

.018 

. 022 

.023 

.018 

* 

. 018 

.009 

.016 

.013 

.018 

.011 

— 

— 

— 

(Transition  downstream) 

F; 

.0  44 

.068 

.204 

. 004 

.0049 

Fc 

. 502 

.680 

.495 

.0385 

. 0201 

ap,/d! 

-.11 

-.22 

-.17 

-.03 

+ .17 

A$.  /f' 

.66 

.64 

.54 

.68 

.61 

-.16 

-.34 

-.31 

-.04 

+ .25 

n 

.30 

.36 

.35 

.021 

.0110 

P(fv)«-at  rj 

-.15 

.-A-00 

. -.  32 

.-.04 

.+.50 

Primes  denote  position  downstream  of  shower. 
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3.  53  Example  of  calculation  of  interfacial  depth  in  a 
deep  estuary  dominated  by  entrainment. 

The  salinity  of  the  ocean  off  Alberni  Inlet  is  about 
32  °/oo.  For  purposes  of  the  computation  we  choose  the  salin- 
ity of  the  upper  layer  as  30.6  °/oo  at  the  mouth.  The  channel 
is  of  more  or  less  uniform  width,  about  1500  m.  The  river 
discharge  varies,  but  we  will  compute  a curve  for  a discharge 
of  60  m-^sec"^.  This  corresponds  to  a fresh  water  discharge 
per  unit  width  - 400  cm2/sec.  The  initial  density  dif- 
ference is  taken  as  To  = 0.02^.  At  the  mouthy  - lO^cm^/sec 
and  = 0.001.  The  depth  at  the  mouth  is  therefore  given 

by  the  critical  interfacial  Froude  number  being  unity 

there 


Thus  !P|  = 4.65xl02cm. 

We  may  now  compute  the  interfacial  depths  at 
various  points  upstream  by  means  of  equation  (3.51.10) 

v VL  =c_ 


where  the  cons  tan 


(2r  is  determined 


at 


the  mouth: 


C. 
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A table  of  such  computations  for  a few  selected 
values  of  "A  is  given  in  Table  3.53*1.  Depths  of  the  top 
.layer  from  Alberni  Inlet  are  given  as  well. 

By  differentiation  of  the  equation  (3.31.10)  it 
is  found  that  the  maximum  depth  of  the  interface  is: 

* '*26  ^ 

and  that  this  occurs  at  the  station  where 

y‘ 

In  the  example  we  are  here  computing 

D,  („»)  - «*  - ao°'3 
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Table  3.53.1 


Computation  for 
in  Alberni 

layer  depth 
Inlet 

Actual  layer 
in  Alberni 

depth 

Inlet 

< 

0^  ( n.  ) 

S °/co 

S °/oo 

, 

D^(m.  ) 

Stat, 

.001 

4.65 

31.0 

31 

4.5 

H 

.0015 

5.6 

30.0 

29.5 

2.0 

G 

.002 

5.2 

29.5 

29.0 

2.0 

F 

.004 

3.9 

27.0 

28 

1.5 

E4 

.012 

2.3 

16. 6 

26 

4.0 

D 4 

.020 

1.8 

6.4 

19 

l 

2.0 

C 

. 025 

. 5 

0.0 

11 

2.2 

B 

5 

2.0 

A 

^Station  is  in  a widening  of  the  channel. 
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3. 6 Salt  wedge 

An  interesting  special  case  of  gradually  varied 
flew  is  the  salt  wedge  which  occurs  in  channels  where  the 
depth  is  insufficient  to  permit  the  salt  water  to  extend 
throughout  the  whole  length  of  the  channel,  30  that  the  down- 
stream part  of  the  channel  has  two  layers  of  fluid  in  it,  and 
the  upstream  part  of  the  channel  has  only  a fresh  water  layer. 
The  wedge  appears  to  be  a phenomenon  in  which  friction  is  more 
important  than  entrainment.  Moreover,  it  is  important  to  dis- 
tinguish between  steady  salt  wedges  and  the  type  of  cold 
front  discussed  by  von  Karman  (Non-linear  engineering  prob- 
lems. Dull.  Am.  Math.  Soc. , 46,  8,  pp.  615-683).  in  von 
Karman's  wedge  vertical  velocities  are  important,  the  pressures 
are  not  hydrostatic,  and  the  slope  of  the  wedge  surface  is 
nearly  1:1.  The  slopes  of  stationary  wedges  occurring  in  es- 
tuaries are  likely  to  be  more  of  the  order  1:40  or  1:100. 

Tiie  hypothesis  which  is  used  in  this  section  to  con- 
struct a theoretical  model  of  the  3alt  wedge  phenomenon  is  that 
the  upper  layer  is  turbulent,  whereas  the  wedge  itself  is  lam- 
inar. '.Ve  shall  use  equation  0.1-16)  which  is  the  equation  of 
motion  of  the  top  layer  vertically  integrated  in  which  there 
is  no  entrainment,  no  free  surface  stress,  and  no  change  in 
section  of  the  channel. 
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The  turbulent  stress  on  the  interface  is  taken  to  be 

(2, 

The  fluid  in  the  wedge  itself  is  presumed  to  be  governed 
by  the  Navier-Stokes  equations,  but  the  velocities  are  so 
small  tuat  the  inertial  terms  are  neglected.  (It  will  be 
noticed  that  the  inertial  terms  are  not  neglected  in  the 
integrated  form  of  the  equation  for  the  top  layer).  The 
equation  for  the  lower  layer,  then,  is  of  the  following  form: 

bV  a 


If  it  is  assumed  that  the  vertical  velocities  are  small 
and  that  therefore  the  pressure  is  hydrostatic,  equation 
(3)  may  be  written  in  the  following  form: 


The  stress  exerted  by  the  upper  fluid  on  the  lower  fluid 
and  given  by  equation  (2)  must  be  balanced  by  a viscous 
shearing  stress  due  to  the  lower  fluid  acting  on  the  upper. 


* 


(5) 
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Me  may  assume  that  the  velocity  used  for  (2)  in  the 
lower  layer  is  given  by  the  following  form: 

v r Cl  +-  -t  G-  ^ 

in  Wi\ich  case  equation  (4)  becomes  the  following: 


and  equation  (5)  becomes  simply 

•Zh.  ^ + 2 fl.  it 

r 


If  the  velocity  of  *,he  lower  layer  vanishes 
at  the  bottom,  CL  vanishes. 

Since  we  assume  no  transfer  of  salt  water 
across  the  interface  there  must  be  no  net  flow  of  salt 
water  across  any  section  of  the  wedge. 

*U.L  pit  ~ o 


(6) 


(V) 


(8) 


Evaluation  of  this  integral  results  in  the 
following  relation: 


CL 


o 


z. 

3 


(10) 
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From  equations  (8)  and  (10)  we  can  solve  for  the  con- 
stant ^ ana  subs ti tuting  this  into  equation  (7)  we 
obtain  the  following  form  of  the  equation  for  the  lower 
layer: 


Llaking  use  of  the  fact  that 


Q - v.d, 


and  equations  (2),  (1),  and  (11),  we  can  eliminate  all 
the  variables  except  J"  and  X so  that  the  enuations 
are  now  reduced  to  two  equations  of  the  following  form: 


t~ 


where  for  convenience 

F = 


(id 


(12) 


(13) 


(U) 


(15) 


(16) 
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We  place  the  origin  of  the  coordinate  system  at  the  head 
of  the  wedge  and  introduce  the  quantity  which  is 

given  by  the  equation 

5,  = D.  of.  A = o 


/ i r-»  \ 

l x ' ) 


We  then  introduce  a change  in  variable 


^ t _ ^ X-i 

a x ~ ° x 


^ , D.  & 


> 2 X 


(18) 


and  rewrite  equations  ( i/+ ) and  (15)  in  the  forms 


/ J_  _ A i ^ 

\ r J d x 3 x ~ 


(19) 


L Xu 
c<  A,  ■ 


< */3  fe)  ■ i 


00 A .3  £ 

D. 


(20) 


Eliminating  the  quantity 
tions,  we  obtain 


3* 


from  the  above  equa- 


( t - A) 


£X. 
<2  x 


c<  A A i_  A , *X-v\ 

D.'l  ^ Ax.  / (21 


.1) 


where  £=  % , the  interfacial  Froude  number.  it 

should  be  noted  that  and  f?  are  not  constants,  but 
arc  functions  of  ( 5,  — ). 
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In  order  to  simplify  equation  (21)  so  that  it 
m iuy  be  conveniently  integrated,  we  make  the  following  sub- 
stitutions which  have  been  substantially  verified  in  the 


laboratory: 


F « I 


A,  = 


/ 


(22) 


liquation  (21)  now 

h)-_  _ 

3 X 


takes  the  form 


(23) 


where  A,5 


the  interfacial  Proude  number 


at  / — O . As  is  a function  of  X only,  equation 

(22)  may  be  written  in  integral  form: 


(24) 


Upon  carrying  out  the  prescribed  integration 
and  collecting  terms,  we  obtain  as  the  final  equation 
which  describes  the  profile  of  the  salt  water  wedge 


Ao(^A  = hr  + A A-.  * (8+**'X/*-i-  (25) 

2 Bo  4 z 
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This  equation  involves  no  other  assumptions  than  equations 
(2)  and  (22  j,  and  will  therefore  be  considered  exact.  It 
is  used  to  determine  the  interfacial  coefficient  of  fric- 


Unf ortunately , equation  (25)  is  very  inconvenient 
with  which  to  work.  We  may  simplify  the  equation  in  trie 
following  manner.  If  (^lo  is  of  the  oruer  1/10  it  may  be 
considered  small  compared  to  the  number  8.  Expanding  the 
natural  logarithm,  in  series  and  collecting  terms  we  obtain 


The  validity  of  this  approximate  form  of  equation  (25)  may 
be  checked  by  comparison  with  the  experimental  data.  jk.  has 
been  substituted  for  J?  in  event  of  a shift  along  one  of 
the  axes  due  to  the  approximation. 

The  interfacial  Froude  number,  /f  , increases 
with  increasing  X • As  ff  cannot  exceed  unity  (see  Sec- 
tion 2.21)  we  may  say 


F £ f=  ± l 

lo  L 


(27) 


The  condition  that ff ~ 1 , therefore  defines  tne  maximum 
tnickness  of  the  wedge,  and  permits  a determination  of  the 
maximum  length  of  tne  wedge.  The  following  quantities  are 
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now  defined 

~ I j X “ X**  } 


(26) 


We  may  determine  the  maximum  value  of  from  the 

eoua  ti  un 


H=i  = 


3 


J/3D/0-O5 


and  therefore 


'6 


tc? 


(29) 


Using  equations  (29),  (23)  and  (13),  we  may  determine 
the  slope  of  the  interface  at  X = Xw, 


fa  oL  f f fu  \ 

/3  * C i - FJn) 


(30) 
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3.61  iSxperimental  studies  of  salt  wedges 

The  experimental  salt  wedges  in  the  laboratory- 
exhibit  a Torn  similar  to  equation  (3.6.25),  as  shown  in 
Figure  3.61.1.  The  solid  lines  indicate  the  wedge  profiles 
determined  by  the  above  equation.  By  fitting  the  experi- 
mental points  to  these  curves,  the  interfacial  coefficient 
of  friction  is  found  to  equal 

Jk-O.  OO30  11 


The  theoretical  curves  have  been  terminated  at  the  value 
of  defined  by  equation  (3.6.29).  None  of  the  experi- 

mental data  exceeded  this  limit. 

In  order  to  fit  the  experimental  data  to  equation 
(3.6.26,  the  approximate  equation,  it  is  necessary  to  make 


- 0.00  s i 


(2) 


Figure  3.61.2  shows  these  results,  the  dashed  line  indica- 
ting the  curve  given  by  equation  (3.6.26).  As  the  approxi- 
mation involved  the  omission  of  a number  of  terns  multiplied 
by  o<~  Ho  it  is  not  believed  that  Jz  bears  any  significant 
relation  to  . 

The  solid  line  in  Figure  3.61.2  indicates  the 
following  empirical  relation  found  prior  to  the  foregoing 
theory 


ki  = <?•  is-(  Ft0) 


,5V 


(3) 
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A further  check  on  the  theoretical  analysis 
would  be  that  the  velocity  in  the  wedge  be  approximate- 


ly correct.  The  maximum  upstream  velocity  occurs  at 


Vi 


WEDGE  PROFILES  A-SERIES  WEDGE  PROFILES  B-SERIES 


u *-  -i  4 


• 8=5* 

• till 


I I I I I 


► ~ O > 


I 


I 

I 


I 


i 

i 
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I ♦ M N ✓ ~ 


Us 


till  li 
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l 
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3.7  'i'he  very  viscous  estuary 
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Certain  estuaries,  though  stratified,  are  so 
turbulent  that  the  deep  layer  is  diluted  v;itn  fresh  sur- 
face water.  These  estuaries  (Type  2)  differ  from  the 
deep-fjord  estuary  (Type  3),  and  river  outlet  (Type  4) 
in  which  the  bottom  iayer  is  undiluted. 

Continuity  principles  show  that  whereas  en- 
trainment can  occur  without  turbulent  flux  across  the 
interface,  turbulent  flux  must  be  accompanied  by  entrain- 
ment. This  and  other  natters  are  discussed  in  Section 
4.6. 

In  this  chapter  the  interesting  feature  of  the 
very  viscous  estuary  is  that  inertia  terms  nay  not  be  im- 
portant as  compared  to  friction  and  pressure  terms.  An 
elementary  theory  of  such  a case  moy  be  developed  as 


f ollows : 


TVio  ^ trno  'nl*  /'ul  nnipfT  nno  ( 1 1 A I o nH  f *5  1 1 7 

a **wk«**X  v u Uw*  v A w**w  ^ »»u  \ y • a.  • X < 


may  be  simplified  to  the  following  forms  if  inertia  terms 
are  neglected. 


(1) 


d y 


* K 


(1 
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where  the  frictional  stress  has  been  written  in  the  form 
~ K(^>  . Jt  is  convenient  to  write  the 

equati  ons  v;i  th  ^ as  independent  variable,  and  let 

V K j(cip  fa)  *hele 

<H'  ' ” ■ * /' 

P 


(2*  ) 


(2’) 


Addition  of  these  two  leads  to  the  simple  re- 


sult that 


sL\i!££*L  +(p-r‘)l£  1 
*i>l  ~L  J 


r O 


(3) 


This  equation.  simply  states  that  in  the  absence  of 
bottom  friction  the  vertically  integrated  pressure  force 
does  not  vary  with 

Now  in  a Type  2 estuary,  with  reasonably  large 
value  of  V , the  term  is  small  compared  to 

the  others,  and  in  many  cases  so  is  ^ ZtV*  it, 
that  the  principal  balance  is  in  the  form 


so 


m.  / 


- O 


U) 


Section  3.7  (Per©  3) 


or  auite  simply 

et  I Q),  -id  l) 

*f'  " ~T‘»F  '"r" 


If  v;c  suppose  that  tiie  density  is  given  by  a simple  law 

- /,  ) 

f 


^ jr»  u.  — 

W X VAX'  ^0 


J 


tnea 


C>j^ 

7- 


- ^ 


V-  " 


From  equation  (4.6.3)  we  can  now  write 


is, 

7- 


s o( 


(6) 


(7) 


New  from  the  dynamical  equation  (2) 


If  we  eliminate 


(g) 


between  ecuatiuns 


(7)  and  (c)  we  obtain  an  equation  giving  £ t in  terms  of 
the  observed  salinity  distribution  (in  a Type  2 estuary) 


i 


‘ ... 

* q x.  * 


19) 
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— This  equation,  despite  its  restrictions  and  its 
approximate  nature,  is  very  remarkable.  If  we  consider 
any  Type  2 estuary,  v;e  can  deduce  the  discharge  ^ , or 
the  non-tidal  velocities,  *'K1  , from  the  salinity  distribu- 
tion, even  though  we  do  not  know  the  friction  explicitly’. 
The  alternate  form  of  (9)  is  as  follows: 


An  interesting  corollary  of  the  equation! 9)  is 


or,  as  in  most  Type  2 estuaries,  if 


(10’  ) 
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3. 71  Sxaaplfc.  Dynamics  of  a Type  2 estuary. 

Pritchard  (1951)  has  discussed  in  detail  the  nu- 
merical magnitude  of  various  terms  in  the  dynamical  equa- 
tion for  the  James  River  estuary  (for  the  discussion  of  the 
salt  transfer  equation  of  use  Janes,  see  Section  4 . 61 ) . i'he 
dynamical  equation  for  the  x- component  is  averaged  over 
time,  a steady  state  is  supposed,  ana  the  following  1‘orrc 
obtained : 


*7  ht+Zrht  = 
>X  7ft 


3£-  * f 3* 


a-x 


£ir 


(1) 


By  analogy  to  the  fact  that  he  found  *7c's'  negligible, 
(Section  4.61)  Pritchard  infers  that  [*’"')  is  also 
negligible.  Hence,  tne  vertical  eddy  flux  of  momentum 
may  be  solved  for  in  the  form 

r 


*1L  w ~ 


I 

uT 


ur 


'tty 


L 


(2) 


i 

Unless  the  elevation  of  the  free  surface  is  known, 

J 0 X 

is  known  only  with  respect  to  a constant  fixed  at  some 
level  . However,  Pritchard  assumes  that  ***** 1 van- 
ishes at  both,  surface  and  bottom  (the  latter  assumption 
being  by  far  the  weakest)  and  evaluates  u1**1  numerically 
from,  equation.  (2).  The  results  for  a single  station  are 
jdnttfld  is  figure  3-71.L, 


Figure  3.71.1  The  vertical  eddy  flux  of 
horizontal  momentum  as  a function  of  depth, 
at  a sample  station  in  the  James  River 
estuary. 
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4. 3 Mixing  processes  in  estuaries 

Very  little  is  known  for  certain  about  the  mixing 
-processes  which  actually  occur  in  estuaries,  so  that  most  of 
the  work  done  to  date  has  been  hypothetical.  These  hypothet- 
ical studies  are  essentially  an  attempt  to  explain  the  obser- 
ved distribution  in  the  estuary  of  a property  like  salt  on 
the  basis  of  assumptions  about  the  mixing.  The  methods  of 
the  tidal -prism  (4,2),  of  mixing  in  segments  (4,3),  and  of 
arbitrarily  defined  eddy  diffusivity  (4.4),  are  examples  of 
this  hypothetical  approach. 

From  the  practical  point  of  view  of  computation 
of  pollution,  numerical  processes  seem  more  likely  to  be  of 
use  than  idealized  hypothetical  mixing  theories.  Section 
4.46  is  an  example  of  the  numerical  process  applied  to  an 
unstratified  estuary. 

An  interesting  feature  of  mixing  in  stratified 
estuaries  is  the  condition  of  "overmixing"  which  is  dis- 
cussed in  Sections  4.51-3,  This  condition  may  turn  out  to 
be  a very  useful  one  in  the  study  of  estuaries  because  it 


does  not  depend  upon  the  detailed  nature  of  the  mixing  pro- 
cess itself. 
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4. 2 The  tidal  prism 

The  first  rough  approximate  method  of  determining 
the  salinity  of  a tidal  estuary  from  a knowledge  of  the  tides 
and  river  flow  is  the  "tidal  prism  method1*,  long-  in  cu»  by 
engineers  (Lletcalf  and  J£ddy , 1935)  in  harbor  studies. 

Denote  the  hush-tide  volume  nf  the  imhmrr  hT 
the  low  tide  volume  by  \JL  , the  volume  of  river  flow  per 
tidal  cycle  by  ♦ The  tidal  prism  is  defined  by 

P - V'h  - VL 

If  the  river  water  is  fresh,  and  the  salinity  of  the 
ocean  is  O , the  salinity  of  the  harbor  is  obtained  by  the 
assumption  that  on  each  tidal  cycle  a volume  f*  of  water  of 
salinity  <J~ ' and  a volume  "R  of  fresh  water,  mix  thoroughly  in 
the  harbor,  and  that  a volume  p + R of  the  mixture  is  ex- 
pelled. The  salinity  of  the  estuary  is  therefore  P c T~  /(p^) 
in  the  steady  state. 
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U.  3 Refinements  of  the  tidal  prism  method 

Ketcnum  (1951)  has  suggestec.  that  the  primary  weak- 
ness of  the  tidal  prism  method  is  that  it  assumes  complete 
mixing  over  the  entire  estuary  during  each  tidal  cycle.  He 
indicated  that  an  estuary  should  be  divided  into  a number  of 
segments  of  a length  approximately  the  local  displacement  of 
the  tide  in  each  of  which  mixing  is  complete.  Before  discuss- 
ing Ke  helium's  theory  of  the  exchange  ratio  we  will  first  dis- 
cuss  various  arbitrary  ways  of  dividing  the  estuary  into  seg- 
ments, some  of  which  may  be  useful  in  cases  where  estuaries 
contain  several  isolated  basins;  and  show  how  dependent  the 
results  may  be  on  the  exact  nature  of  the  mixing  process. 

Suppose  that  the  estuary  be  divided  into  arbitrary 
segments,  in  each  of  which  mixing  is  complete.  The  low-tide 
volumes  of  each  of  these  segments  are  denoted  by  , the 

high-tide  volume  by 

The  segment  ^1  = 0 is  defined  to  be  that  where  the 
quanti ty  ?„  -K*  being  the  river  discharge  per  tidal  cy- 
cle. 

If  ~fn  is  the  fractional  concentration  of  river  water 
in  the  nth  segment  at  high  tide,  then  the  fractional  concentra- 
tion of  sea  water  in  each  segment  is  /“•/»  at  high  tide.  Let 
be  the  fractional  concentration  of  fresh  water  in  the 
nth  segment  at  low  tide. 
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The  rise  of  level  in  segment  H -0  during  flood 
tiae  may  be  assumed  to  be  due  entirely  to  river  flow,  thus 
^ | . Since  there  is  no  flux  upstream  across  the  seaward 

boundary  of  the  Oth  segment,  al30  ^ * / • 

Yfe  now  evaluate  the  flux  of  volume  and  fresn  water 
across  both  the  landward  and  seaward  boundaries  of  the  nth 
segment : 


Ebb  - Landward  boundary 
Seaward  boundary 

Flood  - Landward  boundary 
Seaward  boundary 


Volume  Flux 

On, 

Un-t-l 


Fresh  Flux 

On  —| 

Oflt  i 'f'n 


_(u„- 1?)  ->»- 

- (l»n«  -R)  - -R)?'n+' 


The  freshness  of  the  water  flowing  in  the  flood  tide  is 

/ 


given  by  f 
Let 

= Z 


c * O 


At  nigh  and  low 


High 


Low 


Volume  of 
nth  segment 

Po  + vn 


V 


Total  Freeh 
water  in  nth 

^egment 


(Pn-*V*)K 


The  first  relation  to  be  satisfied  is  that  the 
flux  of  volume  across  each  boundary  be  ^ each  tidal  cycle. 
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This  is  clearly  so  since  at  landward  boundary: 

On  - -R)  a'R. 

and  at  seaward  boundary: 

Un+,  ~ ^ 

The  next  relation  is  that  the  flux  of  fresh  water 
volume  each  complete  tidal  cycle  by  . Thus  at  the  land- 

ward boundary 

Un4-|  ~ (Uh-K)}.!*.  = "R  (i) 

ana  at  the  seaward  boundary 

or*,  \ (1' 

> 

and  ^/nare  in  general  different,  and  equal  only  when 
mixing  is  complete  on  low  tide.  As  a matter  of  fact,  the 
second  equation  follows  from  the  first  by  induction. 

By  conservation  of  fresh  water  the  change  in  fresh 
water  content  during  the  ebb  is: 

~ & 0*  * = L>n  ^n~<  ~ L>n*1 


(2) 
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and  during  the  flood  is 


/„  (P«  +v«)  -pn  t • - (»*-R) 


A number  of  different  results  may  now  be  obtained, 
depending  upon  the  precise  nature  of  the  mixing  process 


supposed. 


Mixing  process  1.  Suppose  no  mixing  occurs  at 


low  tide  so  that  on  the  flood  the  water  flowing  across  each 
boundary  is  unchanged  in  properties.  Formally  this  is  equiv- 
alent to  supposing  that 

fjlt!  * ™ 13) 


Substitution  of  this  into  equation  shows  that  for  this  pro- 

A 

oo  oo  f r>  1 «r  nr»oo^  p of  of  o f O ^ o r»no  ^ r *•»  4 />  V"s  V I 

VW^W  WiiU  uo  ▲ U UUUUJ  ^ VU  WV  4.W  VUV  a.11  MUJ.  vU  | fj  * 

Mixing  process  2.  Suppose  mixing  occurs  at  low 
tide  and  the  segments  are  quite  large  so  that  in  the  limit 


If  there  is  a certain  value  of  /Tt  , say  /7t  * which 


lies  in  the  ocean  then  it  is  certain  that  there  --- 

The  values  of  at  segments  further  upstream  may  be  cal 


“ 'f’l  ~ 


culated  by  equation  (1). 


-V 
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! 


T ryi  - l 

-|  *A*i  - 1. 

e*/c . 


U) 


There  is  no  difficulty  in  satisfying  the  upper  boundary 
between  segments  0 and  1. 

Mixing  process  3.  .Suppose  all  the  water  that 
moves  on  the  flood  is  ocean  water  in  the  form  of  a wedge. 


then  ^ 
by 


and  the  at  the  ntn  segment  is  computed 


u 


m+| 


(5) 


Mixing  process  U.  Suppose  that  mixing  occurs 
at  low  tide  in  each  segment  as  well  as  at  high  tide,  thus 

There  are  two  relations  available:  equations  (1)  and 

(2),  between  which  and  may  be  eliminated,  thus 

yielding  a recurrence  relation  in  £ and  . 


Sec.  k.  31  (Page  1 ) 


U.  31  Exhange  Ratio 

Ketchum  (1951)  has  introduced  a special  type  of 
segmentation  which  seems  to  this  writer  to  be  the  most 
reasonable  "a  priori"  for  unstratified  estuaries.  In  our 
notation,  it  consists  of  defining  ^ in  terms  of  the  up- 
stream segments: 

V„  = u«  + VC  (] 

Moreover,  Ketchun  has  introduced  a quantity  rV\, 
called  the  exchange  ratio,  defined  in  the  following  wa 7: 

r„  = fh /(T%  (; 

and  has  postulated  that  the  total  amount  of  fresh  water 
in  the  nth  segment  at  high  tide  is  given  by  the  follow- 
ing relation: 

= n/r. 

In  our  notation  this  is  tantamount  to  the  following  deter- 


mination of 
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This  result  is  different  from  the  results  of  any 
of  the  four  examples  of  mixing  processes  riven  in  Section 
4.3.  In  order  to  illustrate  these  differences  a very  simple 
hypothetical  example  is  worked  out  here. 

Example.  Consider  an  estuary  which  has  segments 
71  = 0,  l,  2,  3,  but  beyond  these  is  connected  to  the  sea. 
The  segmentation  is  taken  as  being  that  of  equation  (1) 


P* 

V « 

On 


1 

a 


1 

2 

3 

1 


2 

2 

5 

3 


3 4 

3 Ocean 
7 Ocean 
5 8 


Table  4. 31. 1 


The  method  of  hetchum  gives  the  following  values 
of  ~fn  by  equation  (4). 


n 


0 12  3 4 

1 1/2  1/2  l/3  0 


Table  4. 31.2 


Mixing  process  1 leads  t-o  the  result  that  all 
value.0  of  — | 

Mixing  process  2 clearly  is  not  applicable  be- 
cause the  segments  as  given  by  equation  (1)  are  clearly  not 
of  a size  adecuate  for  the  limit  to  be  reached.  Mixing 
process  3 leads  to  the  following  values  of  by  equation 
(5) 


12  3 4 

1/3  1/5  1/8  0 


C 

1 


Table  4.31.3 
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Ilixing  process  4 leads  to  the  following  equation: 


(5 


The  value  of  "/j  must  be  determined  by  the  equation 
4. 3. 1 , where  ^ ^ ® & 


V#  (vh-  V#“^) 
KV* 

■f* 


0 12  3 

0 3 9 15 

-2048 
2 3 5 7 

1 107/135  8/15  1/8 


Table  4.31.4 


By  using  these  different  hypotheses  regarding  the 
mixing  process,  the  salinity  distribution  in  the  sample  es- 
tuary is  quite  different.  Had  we  used  some  method  of  segmen- 
tation other  than  Ketchum's,  it  would  have  been  different  from 
any  of  the  above.  For  example,  the  mean  of  the  estuary 
by  the  tidal  prism  method  is  1/8. 


\x) 
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U.U  A mixing-length  theory  of  tidal  mixing  (Arons  and 

Stommel,  1951) 

Consider  an  estuary  of  uniform  width  W , depth  H , 
length  L_  . The  origin  is  placed  at  the  river  inflow  where 
the  salinity  is  maintained  at  5*0  . The  y -axis  is  dir- 
ected positively  downstream.  At  the  seaward  end  of  the  es- 


n I 

v/  va  x jr  • 9 ' 


~ 4 r-  4 X 

oui  Xlii.  \JJ 


n«^  4 ^ f ' ^ w*  »n  #3 

Xll  OUX  U 


x W x 
v nd  is  w x 


open  ocean  S = CT“  . 

The  river  discnarge  is  ^ leu.  ft/min).  The  mean 
velocity  of  water  in  the  channel  due  to  the  river  is  there- 
fore  a = D 


If  the  length  of  the  channel  is  small  compared  to 
a quarter  tidal  wave  length,  the  tide  will  be  simultaneous 
and  uniform  over  the  entire  channel,  and  we  may  express  the 
height  of  the  tide  as  /*  A.c~“'r  , where  U>  is  the  an- 

gular frequency  of  the  tide. 

The  tidal  current  U is  obtained  from  the  equation 
of  continuity: 


a y/jt  * - H Wa* 


(J  ^ (Jo  slcw-t 

(2  ) 


where 


ue  * 


V»> 
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The  average  tidal  displacement,  S > is  obtained 
by  integration  of  the  tidal  velocity 


$0 


co 


U) 


wnere 


* “ s» 


■v  /H 


i5) 


We  '..ill  consider  the  equation  describing  the  mean 
salinity  distribution: 


bs/dt  * 'U.ts/i*  * *(A  */**)/£* 


(6) 


In  this  equation  S Is  the  time  mean  salinity  at 
any  point  pv  , 'U.  is  the  time  mean  velocity  at  ^ , wnich  we 
may  take  as  that  portion  of  the  flow  due  to  the  river  (that 
is,  SU*-  3i  ),  and  f\  is  eddy  diffusivity  along  the  -axis. 

We  express  in  terns  of  a dimensionless  number  , a char- 
acteristic velocity  which  we  take  at  the  tidal  amplitude  U0  , 
and  a characteristic  length  which  v>e  take  as  2 > the 

total  excursion  of  a particle  due  to  the  tides 


^ a gB  5.U, 


(7) 
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O0  C * ±4  • *c 


This  fore  of  diffusion  equation  regards  the  tides 
as  a turouient  motion  superposed  upon  the  steady  river  flow 
tnrough  the  estuary.  The  simple  assumed  form  of  the  eddy- 
diffusivity  coefficient  is  the  equivalent  of  Ketchum’s 
assumption  of  the  dimensions  of  the  mixing  volume.  It  should 
be  clear  to  the  reader  that  the  simplicity  of  both  of  these 
formulations  results  from  a certain  vagueness  about  the 
physical  process  involved;  the  effects  of  stratification, 
stability,  vertical  mixing,  bottom  roughness,  and  other  in- 
fluences are  not  investigated. 

In  the  steady  state  the  time  derivative  vanishes, 
and  (6)  is  integrable 

as  » A + CL 

(8) 

where  £ is  a constant  of  integration.  At  ^•O  , both 
and  A dt  /d * 6 t because  there  can  be  no  transfer  of 
salt  up  the  river  by  eddy  diffusion;  thus  the  constant  of 
integration  Q **  0 . 

From  (3),  (5),  and  (7)  we  see  that  A may  be  ex- 
pressed as  a function  of 

A - ZB !?•**'/ H'  (9) 

It  is  convenient  to  introduce  a dimensionless  par- 
ameter y 

a c / L- 


(10) 


Sec.  U.U  (Page  U) 


to  express  the  distance  along  the  channel  in  fractions  of 
the  total  length,  and  a dimensionless  parameter 


F*  a.  H x/z  8 C1*  L- 


which  we  may  call  the  "flushing  number". 

Making  these  various  substitutions  in  (8)  the 
following  equation  is  obtained: 

p*s  . "V  ols/el? 

Integration  by  separation  of  variable  yields 
the  following  expression: 

kx,  ^ / , S * f**  , so  that  the  constant  of  integration  C ^ 

is  given  by 

£ 1 m p"  + -^H. 


Therefore,  it  is  possible  to  write  the  ratio  of  mean  sal- 
inity to  the  ocean  salinity  is  exponential  form 


f ( i - i/a; 

e 


The  family  of  curves  on  the  relation  to  s/<r 
is  shown  in  Figure  U.U. I for  the  various  values  of  the 
flushing  number  . 


Empirical  data  for  both  Alberni  Inlet,  Vancouver 
Island,  and  the  Raritan  Riber,  New  Jersey,  are  plotted  on 


Sec . 
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this  figure. 

The  family  of  curves  is  interesting  for  several 
reasons:  (1)  There  is  a toe  to  the  curves  near  ^ sO 

[2)  The  curves  are  very  sensitive  to  P in  region  0.1  ^ F<  3 0. 

(3)  There  is  a point  of  inflection  at  ^ • 

Convenient  alternative  forms  of  the  flusning  number 


I- « 


PH*- 


D HlT 


RbT_ 

vneji1- 


where  is  the  river  discharge,  “f  is  the  tidal  period,  y i: 
the  clean  total  estuary  volume,  and  is  the  river  discharge 
per  tidal  cycle. 

The  curves  presented  here  were  developed  for  a 
very  much  idealized  situation.  For  that  reason  it  is  some- 
what surprising  and  encouraging  to  find  that  empirical  data 
from  actual  surveys  can  he  plotted  on  the  family  with  such 
good  agreement. 

An  attempt  to  calculate  the  proportionality  fac- 
tor 13  from  the  data  v;as  unsuccessful,  the  values  being  of 
an  erde^  of  magnitude  different  for  the  two  cases.  There- 
fore, it  appears  that  although  the  shape  of  the  theoretical 
curves  3s  in  good  agreement-  vit-h  the  observations  at  hand, 
an  a priori  calculation  of  the  flushing  number  is  not  yet 
feasible.  Nevertheless , the  flushing  number  may  be  a con- 
venient concept  to  characterize  estuaries,  just  as  the  fam- 
ily of  curves  themselves  is  a convenient,  semi  -err.pl  ri  cal  ex- 
pression of  the  mean  salinity  distribution. 
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If  the  mouth  of  the  estuary  is  not  taken  as  I 

but  some  point  >,  (o<-  3,  < i)  further  upstream  is 
chosen  as  the  end  point,  then  the  graph  may  still  be 
drawn,  but  the  flushing  number  is  different.  The  sa3 inity 
S,  at  is  given  by 


5-,  /<r~  - <2 


<>~k) 


.f  a new  running  variable 


is  intro- 


duced 


The  result  of  changing  the  location  of  the  "mouth"  is  that 
the  salinity  distribution  is  still  of  the  same  family  of 
curves,  but  with  a different  flushing  number  * F/>,  • 
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4. 41  Mixing  in  estuaries  dominated  by  evaporation  and 
prec ipi tation. 

Consider  an  idealized  estuary  as  described  in  the 
previous  section.  If  (•  is  the  evaporation  from  the  surfac< 
in  ciu/5ec~x,  the  rate  of  change  of  salinity,  5 , will  be 
given  by:  c 


In  the  case  of  precipitation,  £ is  negative. 


A current  is  set  up  to  compensate  for  the  evap- 


orated water 


- £*/H 


The  diffusion  of  salt  is  governed  (as  in  equa- 
tion 4.4.9)  by  a coefficient  of  eddy  diffusion 


A ■ k* 


wnere 


^ " / H L 


The  steady  state  transfer  of  the  salt  transfer 
equation  is  given  therefore  by 


L oL  yu  i + 5=0 


\ 
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This  may  be  simplified  bo  the  form 


^ -r  c*4-  ^2.  + 4-s  = o 

eu.'- 


where 


0L  - * 


rC  hi 


4-  - 


i.  i—t 
A.  I / 


This  is  a form  of  Euler's  equation. 

f- 

By  substitution  S^*  £ 


and  eliminating  the  common 


r 


factor  we  obtain  the  algebraic  equation 

^ ( fc-i)  +•  a. p-  + * * 


the  roots  of  which  are 


r 


~+, 


The  solution  of  the  differential  equation  is  then 

s - G,  * + Ci  * ’ ^ 


The  boundary  conditions  are  that  at^(*  L.  , 

v 

and  also  that  the  net  salt  flux  vanishes,  i.e.  Sr5“  and 

a*  y,-L 


or 


O S * 
GT  r 

-*r-  * 


A s 


0,L  1 * (Ly.!— 


-4- 


C,  L 1 CL  L 


— ih 


Let 


D,  - C.L 


- / 


s *—2. 


The  equations  which  determine  D,  and  J) 


are 


•r-  * D,  t 

^r-  = JD, 


from  wnich  we  obtain 


The  solution  is 


- o , 


a.  (J  ’ 


In  Figure  4.4.2  this  function  is  plotted  against 
for  various  values  of  . The  quantity,  plays  the 

role  of  the  flushing  number  F"  for  estuaries  whose  salinity 
is  governed  by  a balance  between  tidal  flushing  and  evapor- 
ation. 

* 


EH 


ative 


Figure 


that 


4.4.3  shows  the  function  plotted  for  neg- 
is,  for  precipitation. 
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4. 42  Uniform  mixing  along  an  estuary 

If  liie  mixing  in  an  estuary  is  not  caused  by  the 
tides,  but  by  the  wind,  for  example,  the  eddy  diffusivity  /A 
is  not  likely  to  be  given  by  a quadratic  law,  but  may  be  con- 
stant. Under  these  conditions  the  salt  transfer  eauation  is: 


el* 

*4- 


. 4 


ds 


the  estuary^  S 


8lS  - A 


The  integral  of  this  equation  is: 

kL 


S 

rr 


- a 


(*-0 


(i) 


The  boundary  conditions  are  that  at  y * £ , the  mouth  o? 

and 

ds 


(2) 


(3) 
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4.43  hiring  in  a Pillsbury- type  estuary 

As  shown  in  Section  6.5  there  is  a particular 
snape  estuary  in  which  the  depth  is  constant,  and  the 
breadth  is  cf  an  exponential  shape  in  waich  the  tidal  vel- 
ocity is  independent  of  X , tne  position  up  and  downstream, 
many  natural  estuaries  approximate  this  shape.  It  is  of  in- 
terest to  extend  the  theory  of  tidal  mixing  curves  of  Sec- 
tion 4.4,  which  apply  to  a ciianr.el  of  uniform  cross-section, 
to  the  special  case  of  the  Pillsbury- type  estuary  channel. 

First  of  all,  the  eddy  coefficient  /{  is  constant 
in  a Pillsbury-type  basin  because  both  'ZC  and  5 are  both 
independent  of  'X 


Thus  the  relation  given  in  equation  (4.4.7)  is  changed  to 
the  following  constant  form: 

A = if  B ntj-  / co  , 


The  total  flux  cf  salt  across  a transverse  section  is  given 


(2) 


where 


4-  - £ 


( see  Section  6. 5 ) 
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and  the  other  quantities  are  defined  in  Section  4,4.  The 
local  uon-tiaal  velocity  d.  is  given  by  CL  - 
where  o & is  the  river  discharge. 

Equation  (2)  thus  becomes 


AM.  % 


This  integrates  directly  to 

* 

S 


= e 


40<m4^j  $ 


(i- 


ajr-^o  y 


(3) 


Tnis  nay  be  written  in  a form  exactly  the  same  as 
equation  (4.4.12) 


*=-  e 


R <" 1 -p) 


but  in  this  case  ^ ( - 


and 


*0- 


o&-L 


4 4 S uJ  L , H Qua  •€,  fa 4 ; 

where  L~  j is  the  distance  upstream  from  the  mouth  where 
4-  ~ 4-*  /&  or  ^ is  the  distance  upstream  where 
Thus  we  may  use  the  same  family  of  curves  as  shown  in  Fig- 
ure 4.4.1  for  the  Pillsbury-shaped  channel,  only  F7  and 
^ are  redefined. 

For  a quick  evaluation  of  the  Fj  it  is  convenient 
to  locate  the  value  of  /4r0  where  s/<r-  has  some  given 
value,  say,  0.5.  luble  4.43.1  gives  the  value  of  F,  for 
various  values  of  /i,  at  S'  /y" 


* a.r 


Section  4.43  (Pa^e  3) 


TABLE  4.43.1 


X, 

FT 

« 

0.9 

6.  3 

0.  8 

2.  8 

0.  7 

1.0 

C.  6 

1.04 

0.5 

.69 

0.4 

.46 

0.3 

.21 

0.2 

.17 

0.1 

. 08 

0.  05 

. 037 

0.  02 

.014 

0.  01 

. 007 
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uu  nxnmplc  oi  .horizontal  mixing  theory : me  Severn 

The  data  available  for  the  Severn  Estuary  (7.1J) 
in  adequate  to  test  the  various  theories  proposed  in  the 
previous  sections.  First  of  all,  it  is  clear  that  the  Sev- 
ern is  an  unstratified  Type  1 estuary,  to  which  the  ideas 
of  horizontal  mixing  ought  to  apply,  if  they  are  correct. 

We  will  ;omp”te  horizontal  salinity  distributions  on  the 
basis  of  Ketchum's  Exchange  Ratio,  and  on  the  basis  of  the 
Arcns-Stommel  theory.  It  will  be  shown  that  neither  of  these 
methods  work  for  the  Severn.  The  inference  appears  to  be 
that  the  mixing  length  involved  is  not  even  remotely  similar 
to  the  tidal  displacement  in  the  Severn,  but  is  more  nearly 
the  depth. 

The  Severn  by  method  of  Section  4.31 

The  volume  of  water  in  the  Severn  Estuary  was  com- 
puted from  the  data  given  by  Gibson  (1933).  Figure  4.44.1 
snows  the  accumulated  volume  of  water  from  Gloucester  to  sea- 
ward at  high  and  low  water  during  Spring  tides.  For  torch, 
1940,  the  average  river  discharge  (Figure  7.11.2)  is  2,600 
cusecs,  which  is  equivalent  to  1. 2xl0^f t^/ti dal  cycle.  The 
segmentation  of  the  estuary  is  started  by  recalling  the  defi- 
nition of  the  /JV-O  segment , "R,  = R , and,  using  equation 
(4.31.1)  the  subsequent  and  V,  volumes  may  be  readily  de- 
termined. The  following  table  summarizes  the  calculation: 


- r - • ‘ ' • T *>orr  “>  ! 


TABLE 

U - 4 4 ■ i 

0 

I 

2 3 

Ocean 

xl0"J  - 

1.2 

89 

910 

w 

xJG"^  = 

2.  2 

3.3 

92  1000 

4-t' 

3. 

0.01 

0. 001  - 

0 

0 

0 

11 

0 

32 

32 

32 

0 

0 

0 

0 

8 

16 

32 

From  the  results  through  Segment  2 further  calcu- 
lation to  the  ocean  is  not  warranted. 

The  freshness,  -fn  , and  salinity,^,  are  the 
average  values  for  the  individual  segments.  For  compari- 
son, the  salinities  from  Figure  7.11.5  during  winter  spring 
tides  are  given.  It  is  evident  that  the  salinities  computed 
by  the  method  of  Section  4. >1  are  greatly  in  excess  of  those 
observed. 

The  Severn  by  method  of  Section  4,43 

It  is  easy  to  see  that  Figure  7.H  5 may  be 
plotted  on  a flushing  number  graph  (Figure  it.lt.l).  When 

this  is  done  it  is  found  that  3 is  of  the  order  of  raagni- 

_ 3 

tude  cf  10  , for  both  summer  and  winter.  The  tidal  excur- 


sion is  of  the  order  of  150,000  feet  in  the  Severn.  The 
mixing  length,  therefore,  is  no  more  than  150  fact.  As  a 
result  the  hypothesis  upon  which  Sections  U U to  is 


TOTAL  VOLUME  (FT)3xlO“ 
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FIG.  4.44.1 
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based  is  3iraply  not  satisfied.  The  mixing  length  appears 
to  be  more  certainly  related  to  the  depth.  Upon  reflection 
the  reader  will  probably  agree  that  this  seems  more  reason- 
able anyway.  Instead  of  the  mixing  being  done  by  huge  hori- 
zontal eddies,  several  miles  in  length,  the  mixing  is  done 
by  small  "boils"  or  eddies  resulting  from  the  shearing  flow 
over  rne  bottom.  In  order  to  demonstrate  this  fact  more 
clearly  we  have  considered  the  Severn  in  more  detail. 

The  steady  state  distribution  of  salinity  is  a bal- 
ance of  diffusion  and  advection: 


as  - A — 


where  the  / -axis  is  directed  downstream,  S is  salinity,  ^ 
is  the  mean  river  velocity,  and  A the  coefficient  of  ciffu- 
sivity.  The  quantity  any  X may  be  computed  from  the 

river  discharge  Q divided  by  the  width  and  depth  ^ of 
the  estuary: 


- a luM. 


The  diffusivity  may  be  written  as 

A = 13  lcLnr 

where  or  is  the  amplitude  of  the  tidal  velocity  and  is 
the  depth. 

In  this  case 

1 - ( Qs  'j  j ( J 
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R 

We  computed  the  value  of  O at  ceveral  points  in 
the  Severn  from  this  formula. 

TA3LE  4.44.2 


« 

/eston  Port 1 ahead 

A u s t 

Sh?rp.n.?^Cf 

A t*1  r\  /«*  V-» . 

* A-i.  A JL 

(p*  (feet)xlO' 

46.  0 

26.  0 

6.9 

5.2 

2.  1 

qL.  (feet) 

70 

60 

50 

20 

15 

f\j*  (ft/sec) 

P c 
w ° / 

2.  5 

8.  5 

8.  5 

8.0 

Winter 

S'  (°/oo) 

23 

16 

8 

6 

4 

oL%  jd>fL  (°/ooxioVft) 

0.6 

0.8 

0.  8 

1.0 

1.2 

Summer 

s (°/co) 

28 

27 

25 

20 

18 

/ooxloVft) 

0.2 

0.2 

0.  2 

0.8 

0.6 

“3  f winter 

0.6 

0.7 

1.  5 

9.0 

5.0 

3 ^ summer 

0.  3 

0.7 

3.0 

5.0 

6. 0 

Winter  Q ftVsec 

2600 

Summer  0 ft-^/sec 

360 

The  Raritan  River 

The  Raritan 

River  is 

the  chief  example  given  by 

Ketcnun  (1950)  to  illustrate  the  hypothesis  of  the  exchange 
rafio.  It  is  interesting  to  determine  the  mixing  length  at 
several  stations  similar  to  those  discussed  by  him 
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Station  (1/  3 miles  upstream  of  South  Amboy 
Station  (2)  South  Amboy 

Station  (3)  3 miles  downstream  of  South  Amooy 

The  river  discharge  used  by  Ketchum  is  33xlC^f t^/12.  5 
nours,  or  730  ft^/sec.  The  tidal  velocity  is  about  1.5  knots, 
or  2.6  ft/sec. 


Station  1 

Static  n ? 

Sta t ? on  3 

Depth  (ft.  } 

9 

12 

18 

Width  (ft. ) 

3600 

7800 

24,000 

S °/oo 

10 

25 

26 

JU  /y  °/ oo/rai 

4.4 

0.4 

0.2 

pu/ty  °/oo/ft 

0.7X10-3 

0.7x10"^ 

0. 3xlO"4 

Mixing  length 
in  feet 

128 

1100 

560 

The  mixing  length  is  clearly  much  greater  than  the 
depth,  but  not  as  large  as  the  tidal  excursion. 

The  reader  will  see  that  whereas  the  mixing  length 

i 

in  the  Severn  is  of  the  order  of  magnitude  of  the  depth,  in 
tne  Raritan  it  is  much  larger,  but  not  nearly  as  large  as 
the  tidal  excursion. 
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4.  45  Horizontal  tidal  exchange  through  an  inlet. 

A constricted  inlet  may  be  expected  to  act  as  an 
efficient  tidal  exchanger  because  of  the  tendency  of  the 
flow  into  it  to  be  potential,  whereas  the  flow  out  is  jet- 
like, as  indicated  in  Figure  4. 45.1.  As  a result,  one 

should  expect  that  for 
the  most  part,  the 
water  passing  through 
the  inlet  when  the 
current  reverses  is  not 
the  same  as  that  before 
reversal. 


pic  ture : 


Let  us  consider  the  following  simple  theoretical 
an  inlet  of  width  A (Figure  4.45.2).  For  sim- 

CK. 

plicity  we  assume  that 
the  depth  ^ is  uniform 
both  in  the  inlet  and 
throughout  the  estuary. 
The  discharge  on  the  ebb 
is  from  a semi-circular 
region  of  radius  , and 
amounts  therefore  to 


OC  <»  w 


Figure  4.45.2 


«D  . The  flow  during  flood,  which  occurs  in  a jet- 
like  filament  of  width  CL  , and  length  , amounts  to  — ^ 4£> 
The  total  river  discharge  for  a complete  tidal  cycle  is  , 
and  for  simplicity  is  written  as  -R  ^ arp 


T 

I 
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The  conservation  oi  mass  requires  that 

=ar  (1) 

The  conservation  of  salt  requires  that 

+!^^-cJr)s  =0  (2) 

where  5 is  the  salinity  inside  the  inlet  and  .ns  the 
salinity  of  the  ocean.  The  assumption  involved  is  that 
on  the  ebb  the  portion  of  the  jet  within  the  semi- 

circle ebbs  undiluted,  but  that  the  remainder  of  the  ebb 
is  water  of  salinity  $ . Elimination  of  ^ between 
equation  (1)  and  (2)  results  in  the  following  expression: 


X 
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4.46  Computation  of  pollution  in  a vertically  mixed  estuary 
by  numerical  process . 

As  was  seen  in  Section  4. 44*  there  appears  to  be 
serious  difficulty  in  applying  to  an  estuary  any  of  the  hypo- 
thetical mixing  processes  discussed  in  the  early  sections  of 
this  Chapter.  For  example,  the  salinity  distribution  in  the 
Severn  estuary  did  not  seem  to  fit  properly  unat  computed  by 
the  exchange  ratio.  The  method  of  Section  4.44  is  always 
open  to  the  objection  of  the  restrictive  assumption  about  the 
geometry  of  the  estuary  and  the  undetermined  nature  of  the 
constant  B.  From  a practical  point  of  view  the  proper  procedure 
is  to  use  the  distribution  of  river  water  as  a means  of  discover- 
ing the  magnitude  of  the  turbulent  diffusion  coefficients  at 
various  places  in  the  estuary,  and  to  devise  a method  using 
these  coefficients  which  will  yield  the  dilution  of  pollution- 
at  any  point  in  the  estuary.  An  attempt  at  such  a method  is  de- 
scribed in  the  following.  It  is  important  to  re-emphasize  that 
it  is  intended  to  apply  only  to  vertically  unstratified  estuaries 
in  which  the  mixing  is  due  to  tides. 

The  steady  state  equations 

Let  us  consider  an  estuary  such  as  that  shown  in 
Figure  4.46.1.  The  ^ -axis  is  directed  along  the  axis  of  the 
channel,  the  cross-sectional  area  S(*)  may  vary  with  the  po- 
sition along  the  axis. 

Now  we  will  suppose  that  a pollutant,  which  is  mis- 
cible with  water,  whose  average  (over  a tide)  concentration  c (yj 
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varies  with  , is  in  a steady  state  distribution  in  the 
estuary. 


The  term  steady  state  means  that  the  average  of 
the  concentration  (2,  over  a tidal  cycle  does  not  change  from 
tide  to  tide.  This  will  be  true  if  there  has  been  little 


nhanup  in  "hhp>  T*i 'tr^T*  rl  i srbafjp  diiri  n~  the  ti  ms  i'"'VQlvpa 

O- 0“  . w-1. 
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the  discharge  of  pollutant  into  the  estuary  has  remained 
constant.  If  the  total  river  discharge  into  the  estuary,  at 
some  point  remote  in  the  — J(  direction,  is  Q'.  then  the 
flux  of  pollutant  by  advection  toward  the  sea  is(2c. 

In  addition,  there  is  a turbulent  flux  due  to  the 
tidally  produced  turbulence  in  the  estuary.  This  turbulent 
flux  may  be  written  in  the  customary  fashion;  -5f\ 
where  A(x)  is  a turbulent  eddy  diffusivity,  the  value  of 
which  must  be  determined  before  the  distribution  of  the 
pollutant  can  be  calculated.  Contrary  to  the  methods  of 
hetchum,  and  of  Arons  and  Stommel,  we  will  not  specify  /A 
a priori. 

The  net  seaward  flux  FO*)  of  pollutant  across 
any  section  is  the  sum  of  these  two  fluxes. 


FM 


(1) 


If  the  pollutant  is  conservative  (does  not  decay  with  time), 
the  net  flux  must  be  constant  downstream  of  the 

source  of  pollutant,  and  must  be  zero  upstream  of  it.  Some 


i 


igure  4.46.2 


Section  4.46  (Pace  3) 


) 

1 

» 

important  pollutants,  such  as  bacterial  ana  a conic  pile 

! 

wastes,  are  not  conservative,  and  their  concentrations  de- 
cay even  v.hen  not  subject  to  the  dispersing  influence  of 
tne  estuary. 

1 • 

The  concentration  of  such  a non-conservative 

t 

pollutant  in  an  isolated  container,  nav  decrease  with  time 
according  to  the  exponential  decay  lav/ 

— 1 1 / ~c 

c = c0  e 

where  £ is  the  time  required  for  the  concentration  to 
decay  from  the  concentration  Co  to  . The  time,  ~C  , 

is  slightly  larger  than  the  half-life  of  the  pollutant: 

0.693  T - (half  life) 

The  net  flux,  F(v)  , cannot  be  the  same 
for  all  values  of  X if  the1  pollutant  i3  non-conservative, 
but  must  diminish  in  the  following  manner: 


The  steady  state  equation  (1)  may  be  written  in 
the  following  form: 


at  outfall 
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where  fy/  is  the  total  rate  of  pollutant  supply  at  the 
outfall. 

For  a given  pollutant  and  gjven  estuary  it  is 
evident  that  the  quantities  ^ , and  T.  are  known. 

If  v;e  can  determine  , then  the  distribution  of  c(x) 


K t . 1 ^ V.  ~ 4 1— 1 J 
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Determination  of  the  eddy  dlffualvlty  “ 

The  eddy  diffusivity  A may  be  computed  from 
a knowledge  of  the  concentration  of  any  other  property  in 
the  6stuary.  The  concentration  of  fresh  water ) , i3 
a convenient  conservative  property.  In  some  estuaries 
v/  here  ground -water,  or  precipitation,  or  evaporation,  is 
important,  it  might  be  almost  impossible  to  use  the  con- 
centration of  fresh  water  as  a means  cf  computing  A . 

Let  us  suppose  for  the  immediate  purpose  that  all  the  fresh 
water  comes  from  one  river,  and  that  little  is  added  or  sub- 
tracted by  other  means.  How  to  treat  more  complicated 
cases  will  be  evident  when  the  principles  for  this  simple 
case  are  understood. 

If  the  concentration  of  fresh  wa^er,  , is 
written  for  in  equation  (1)  we  obtain  the  following  form 
of  the  equation,  because  in  this  case  FM-a 


S'  cLl/djC 


(5) 
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Tile  concentration  of  fresh  water,  “f"  , is  dimen- 
sionless  (e.g.  , pure  fresh  water  is  represented  by  1.0;  a 
mixture  of  one  part  fresh  water,  three  parts  sea  water,  is 
given  by  = 0.25,  etc.).  The  quantity ,_ Q , (river  dis- 
charge)* has  dimensions  of  ft^/sed;  the  Cross^-'sectionaT  "area , 
ft^;  and  distance  in  the  X direction,  ft.  Thn  dimensions .... 
of  the  eddy  diffusivity  A are  therefore  ft^/sec. 

-i-.  = The  easiest  way  to  compute  A is  to  segment  the 
estuary  along  its  axis , the- intervals  being  equally  spaced 
0k  -feet  apart  (Figure  4.46.2 JY 

The  average  value  of  is  indicated  at  each 

segment.  The  equation  (5)  may  now  be  written  in  finite  dif- 
ference form  for  numerical  computation: 

_ . _ Q2a.  / / -4^) 

A - = S«(4»- 7- W l6) 

In  generaT,  is  less  than  unit,  and  ~fti  +1  is  less  than 
, so  that  Ah.  is  expected  to"he  positive.  The  values 
of  at  points  far  upstream,  (where  -5  / and  -fn-t  ^ ) 

are  indeterminate  from  the  fresh  water  distribution;,  this  is 
also  true  beyond  the  mouth  of  the  estuary  in  the  ocean. 
Calculation  of  the  concentration  of  pollutant 

Equation  (4)  is  written  in  finite  difference  form 
and  the  various  coefficients  gathered  together 


rh  M-* 


r\ 


+ 


*-  rt-f  | 


r O 


(7) 
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Equation  (7)  must  hold  at  each  of  the  segments, 
with  one  exception.  In  addition,  there  are  certain  addit- 
ional constraints.  At  the  ocean,  where  -f- ^ Q , the  value 
of  C~0  , and  far  upstream,  C 6 . At  the  segment  where 

the  outfall  of  the  pollutant  la  located  one  further  condi- 
tion is  imposed:  the  difference  in  flux  upstream  and  down- 
stream must  be  equal  to  the  total  rate  of  input  of  pollutant 

/ IJ  VJ  o tio  u r\  oo-?  ;»nn  f n V->  r>  00  o **»  o +-  t»»  K ♦-  V>  f*  nnppop 

• H w iuka  j na  v/  o j.  uii'/  o v U ta  is  ••  lix  vs  il  1/  ills  U1  V v 

(outfall)  of  poll',  .nt  is  nearest  as 

In  place  of  equation  (7),  which  holds  at  all 
points  except  at'VU^f,  y;e  have 


(7*  ) 


The  equations  having  been  obtained,  the  solution 


is  easily  carried  through  by  relaxation  methods. 
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Simple  example  .a^tmrw  of  uniform  section 

before  proceeding  with  a numerical  example,  it 
may  be  or  interest  to  consider  a simple  example  to  illus- 
trate certain  fundamental  properties  of  the  problem.  We 
will  choose  an  example  which  admits  an  analytical  solution. 


r\  ^ ^ A dJ  ~ ~ 
w uuo  mu  c;  a. 


<3  o u UQl  j 


UUli Uii 


\jx  uao  oou  uiuii 


which  opens  abruptly  into  the  sea  at  ^ ~ O . If  the  eddy 
diffusivity  ^ is  a constant,  then  the  concentration  of 


fresh  water  has  a very  simple  analytical  form: 

/ ' 


| ^ f-  e 


* 


(8) 


where 


-x 


/ - 


_Q 

As 


The  water  is*" entirely  fresh  f ar  upstream  (-£■=.  / ) , but 

diminishes  to  -^»-c  at  the  mouth,  as  shown  by  the  dashed 
line  in  Figure  4.46,3.  In  this  graph  the  distance  along 


the  estuary  is  plotted  in  dimensionless  units  “X 1 ( X is 

negative  within  the  estuary). 

Now  consider  the  concentration  of  pollutant  due 
to  an  outfall  at  *X  ” L . The  analytical  solution  of  equa- 
tion (4)  is  composed  of  two  separate  parts: 


"lies  upstream  of  the  outfall 
plies  downstream  of  the  outfall 
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The  quantity  £a  is  the  concentration  of  pollution  which 
would  result  at  the  outfall  if  there  were  no  diffusion 


or  decay,  in  other  words,  the  simple  dilution  by  river 
flow  alone. 

Fig.  4.46.3  displays  certain  of  these  solutions. 
The  abscissa  is  expressed  in  a dimensionless  fashion. ^The 
ocean  is  located  at  the  extreme  right-hand  side.  Toward 


Figure  i4.46.3 


the  left  is  the  upstream  direrti oh". “ The  dashed' lia-e  is 

the  ,f,he.~eqnp>fri.aa  Jloru_the--^ficent^ti-©&  -of- -fresh: 

river  v/ater  , (equation  8).  The  solutions  of  the  prob- 
lem of  the.„diatri  button  of.  a conservative  pollutant  intro- 


duced at 


ir*  r -3  -2.,  -f  -O.r;  -<5. 

A — / > — = 


are -shown  as  soiicLblack  curves,  terminating  at-  the  dashed 
curve  at  each  of  the  outfall  lojcat ions ^raspectri vely.  The 
remainder-  of  - the  ; - s olu  t i on  ( doivns  tream  of  the  outfall)  falls 
along  the  dashed  iine0  The  peak  concentration  is  at  the  out- 
fall;'  _ It  is  immediately  clear,  then j_in  this  simple  case  of 
an  estuary  of  uniform  cross -section  and  eday-dif f usivity , 

- that-  the  pollution  extends  far  upstream  of  the  outfall,  no 
matter  what  the  location  of  the  outfall,  but  'that  the~  con- 
centration is  everywhere  reduced  upstream _cf  the  outfall  if 
the  location  of  the  outfall  is  moved  toward  the  sea. 


A seawafcT'drsplaceraent  of  the  outfall  also  re-- 
duces  the  peak  cphcentraticn  at  the  outfall,  but  dovm stream 
the  concentration  of  a conservative  pollutant  is  not  affect 
ted  by  such  a displacement. 

The  dotted  curves  'in  Figure'  7+ 746. ’3  represent 
solutions  for -a  -non-conservative  pollutant  whose'  decay  time" 

r * A S H fed  ) - ::r 

for  outfall  locations  at  ^ = -3  ana  -1„  The  remarkable 

feature  of  these  curves  is  that -the  peak  cone entrati bn/^e yen- 
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.0 


at  the  outfalls,  is  much  reduced 
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teak  concentration  is  as  fellows 


fro  si  what  it  is  in  the  case 
rue  general  equation  for  the 


fne  non-conservative  pollutant  extends  both  up  and  down- 
stream of  the  outfall,  but.  unlixe  the  conservative  pollu- 
tant, can  be  reduced  at  a point  below  the  outfall  by  an. 
upstream  displacement  of  the  outfall. 

Consideration  of  this  simple  case  also  shows 
clearly  the  effect  of  a change  in  the  river  discharge  . 
Increase  of  river  discharge  increases  the  fresh  water  con- 
centration at  every  position,  foreshortens  the  extent  of 
pollution  upstream,  reduces  the  peak  concentration  of  pol- 
lution at  the  outfall,  and  even  reduces  the  concentration 
of  a conservative  pollutant  downstream. 

Numerical  example 

As  an  example  of  computation  of  concentration  of 
pollution  by  numerical  methods,  v particular  distribution 
is  computed  for  the  Severn  estuary,  for  which  there  i3  rood 
salinity  data  available  (Bassinaale,  19-+3).  Stations  are 
chosen  beginning  at  Gloucester  ( * t );  the  sections  are 

spaced  about  two  miles  apart  (^,  = 10,000  feet).  Cross- 
section  areas  are  given  in  Taole  k.u6.1,  in  the  column 
headed  . The  fresh  water  concentration,  from 


actual  summer  time  survey  data,  is  also  tabulated. 
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The  eddy  diffusivity,  , is  computed  at  each 
station  by  equation  (6)  using  a value  of  360 -ft3/  sec 

as  the  fresh  water  discharge  during  the  summer  dry  period. 

It  is  interesting  to  note  that  the  highest  values  of  An. 
occur  at  stations  6,  7 and  8,  which  happen  to  be  the  lo- 
cation of  the  Severn  bore. 

Suppose"  now  that  a pollutant  outfall  is  located 
at  Sharpness  (VVL  = 12} , and  that  this  pollutant  has  a half 
—life  of  4 days  ( "C  = 5xlC)5sec).  The  coefficients 
of  equation  (7)  are  tabulated  in  Table  4.46.1.  An  estimated 
distribution  of  pollutant  is  initially  entered  into  the 
Table,  and  the  sum  (called  the  residual)  is  formed: 

^ <Lk.i  + ■*" 


(ID 


Equation  (7)  requires  that  all  residuals  except 
that  at  the  outfall  should  vanish.  In  fact,  we  see  that 
this  is  far  from  true  for  the  estimated  pollution  distri- 
bution. We  can  adjust  the  value's  of  Crl  by  success tive 
steps  (relaxation  method  of  Southwell,  1945)  because  it  is 
obvious  that  if  we  add  b ^ units  to  any  C ^ not  only  will 
the  residual  be  changed  by  the  addition  of  O.AC-i  but 
also  the  residual  will  be  increased  by*]^^  and 

the  residual  by  T?*-,  ac*  r We  make  these  adjust- 

ments step  by  step,  removing  the  worst  residuals  first,  and 
from  time  to  time  re -computing  the  residuals  from  the  newly 
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( 

adjusted  Cn  by  equation  (11).  When  the  residuals  have 
all  been  reduced  to  a point  where  further  read  jus tments  of  C M 
are  below  the  precision  of  the  reouired  result  (say 

| 

l/o  to  5 yv)  the  relaxation  procedure  may  be  regarded  as  com- 
pleted: namely,  we  have  obtained  a set  of  values  of  ^ 

which  satisfy  equation  (7).  in  the  example  shown,  54  ad-  1 

jactnents  of  C*  were  necessary  in  all.  The  computation 
took  about  three  nours.  We  have  used  arbitrary  units  for  j 

C ,y  and  must  new?  convert  to  units  of  actual  pollution 
concentration.  We  may  suppose  that  the  discharge  of  pol-  i 

lutant  at  station  12  is  10  lbs/sec.  The  concentration  due  1 

to  uilution  by  a river  flow  of  360  ft-Vsec  alone  (if  there 
were  no  diffusion  cr  decay)  would  be,  therefore,  28xl0”^lbs/f . 

The  values  of  <ij)  (=  10  lb3/sec)  is  put  into  equa- 
ti un  ( 7 ’ ) , f/z*-  should  then  be  equal  to  the  residual  . 

I 

If  all  the  concentrations  in  arbitrary  units  are  multiplied 

t 

by  5x10  the  value  of  comes  out  properly.  The  values 

of  the  concentration  of  pollution  thu3  computed  are  shown 
in  the  last  column  of  Table  4.46.1. 
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4.51  The  effect  of  controls  on  the  limitation  of  mixing: 
"overmixing". 

There  is  a curious  way  in  whies  a control  acts  to 

limit  the  possible  amount  of  mixture.  Consider  an  estuary  F" 

opening  into  the  ocean  ^wi  thout  tides^  Q through  a transition 

T.  If  the  entire  system  is  sufficiently  deep,  and  the 

O 


river  R discharge  ( ^ % per  unit  width  of  the  transition) 
is  not  too  great,  there  is  a two  layer  flow  through  the 
transition.  The  depth  of  the  upper  layer  is  determined  by 
the  method  of  Section  2.21.  If  there  is  no  mixing  between 
the  two  fluids,  the  lower  layer  is  stagnant,  and  the  upper 
layer  at  *J*  is  entirely  fresh. 

Suppose  that  some  agency  for  vertical  mixing  of 
the  two  layers  exists  in  the  estuary  and  that  the  amount 
of  mixing  is  progressively  increased.  The  upper  layer  is 
now  somewhat  brackish,  the  discharge  of  both  layers  through 
Y*  is  increased,  and  the  interface  nearer  to  mid-depth.  In- 
creased mixing  in  E decreases  the  salinity  difference  of 
the  two  layers  at  T~  , and  increases  the  discharge;  but  there 
is  a point  beyond  which  increased  mixing  has  no  further  effect 
on  either  the  discharges  through  or  the  salinities  at  T . 


It  is  natural  to  suppose  that  some  estuaries  will 
be  over-mixed  in  this  sense,  and  that  the  salinity  control 
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v/iil  be  exercised  by  the  nature  of  the  transition  rather 
than  by  diffusion,  horizontal,  vertical,  or  otherwise. 

ft  * 

The  criterion  for  overmixing" can  be  obtained  from 
the  equation  for  a stationary  interfacial  wave  (2.21.9'): 


*/•  <*✓ 
a / 


t *±IL 
(1-7)'  7' 


the  equation  of  continuity 


and  the  equation  for  the  conservation  of  salt. 


To  relate  salinity  and  density  we  may  take  a simplified 


equation  of  state 


= / +a.s 


Thus 


^3  - f*~fi  Cs"  a.  ( *■*  j 
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The  equation  ( 1 ) may  now  be  written  in  the  following  form 


or 


J (it)'  + i - 
V y 

I1  ] 


(V) 


3*^/, 

l'v 


st  \ 5 /v? 
sZJ  ^ 


(8) 


It  is  convenient  to  introduce 


= ■t'-y + (i-y3  - -n'l-D^yii-iP'o 


19) 


This  is  an  approximation  of  equation  U). 

The  real  roots  of  f(y = ° between  'y  = o and 
- / cease  to  exist  a3  i) — ■*  / and  this  point  determines 
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is  a double  root  occurs  when 
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If  v.e  multiply  equation  (10)  by  ‘Vj  and  subtract  from 
(9)  we  obtain 


jt  (>--*)  7 


•*.*»*=/ 
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Solving  (11)  for  t)  and  substituting  into 


(9)  yields 


This  can  be  conveniently  rearranged  in  the  form 


(12) 


7 


* , (y  v-  Ci-yJ  / 


■/t  3 


(13) 


— . This  means  that  for 


As  ^ o°,  

weak  river  flows  in  estuaries  with  large  transitions  the 

cvermixing  leads  to  an  interface  at  half  depth.  In  order 

to  investigate  approximately  the  effect  of  ^<00  we 

1 

introduce  £ ^ ^ ^ is  a small  quantity  as  long 

as  is  large.  Substituting  this  into  equation  (13) 
and  retaining  only  the  large  terms  we  obtain 


e. 3 - s - i = ° 


(H) 


For  large  values  of  jL  the  roots  of  this  expression 
are  approximately 

3 


£ - 


- ’{Wi 


(15) 


In  this  way  we  can  evaluate  in  states  of  overmixing 
at  what  values  of  the  quantity  deviates  signif- 

icantly from  l/2.  The  ratio  of  salinities  of  the 
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4.  52  Experimental  study  of  overmixing 

Two  sets  of  runs  were  made  past  a .^transition  to 
test  the  theory  of  overmixing.  The  transition  was  a narrow 
segment  of  width  opening  to  full  width  of  the  channel 
both  upstream  and  downstream.  The  length  of  the  narrow  sec- 
tion was  42  cm  in  series  G-,  and  100  cm  in  series  H.  Sample 
runs  are  summarized  in  Table  4.52.1  and  Figure  4.52.1.  In 
all  runs  the  total  depth  Q^was  31.5  cm,  is  the  river 

discharge  per  unit  width  of  the  transition.  A is  the  den- 
sity- of  rTver  disc har ge_;  ^ the -density  of  the  overmixture-;, 
the  density  of  the  sea-water.  is  defined  as 


A 


-0  «--■  A fr 

- hf  p>.  -/•- 


and  ;T^ "is  computed  fro^the  .theory  of  Section  4.  51.  The 
data- -is  plotted'  id  Figure  y4.:52.T^,’--the--G-irG*les  -are--^  -ru-ns-|— 
the 


s are 


series p- the  square  is  the  only  run 


( Qj-  ~t>  ) where-  £,  is  of  sensible  order  of  magnitude. 
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Overmixing  experimental  runs 
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P 
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FIGURE  4.52.1 
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4 . '5 3 Influence  of  tidal  currents  on  overmixlng 

Pour  tests  were  made  In  which  a tide  was  imposed- 
on  the  overmixed  condition.  Figure  4.53.1  shows  the  result 

-“,-=  7"  ~ “ . • » '■raair.y-.' 

of  these  experiments.  In  test  2,  with  the  two  feute  tidal 
period,  and„test  3,  the  sea  W8ter  of  density  is  swept 
out  of  the  transition  for  a portion  of  the  tidal  cycle  and 
there  is  a reversal  of  flow  in  the  upper  layer.  It  is  evi- 
dent that  serious  disagreement  in  the  of  the  steady 

and  unsteady  state  tests  is  not  experiences  until  test  4.‘ 

In  this  test  the  surface  profile  was  distorted  by  a bore 
and  there  was  no  vertical  stratification. 
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4.  54  Overmixing  in  actual  estuaries 

The  only  type  of  estuaries  to  v*hicn  the  concept 
of  overmixing  can  properly  be  applied  are  those  which  are 
vertically  stratified  and  in  which  there  is  evidence  of  a 
net  flow  (averaged  ever  a tidal  cycle/  downstream  in  the 
upper  layer,  and  a net  upstream  flov:  in  the  bottom  layer, 
it  is  important  not  to  atter.pt  to  apply  this  hydraulic  con- 
dition to  vertically  mixed  estuaries,  such  as  the  Raritan 
River,  Delaware  Bay,  Bristol  Channel-Severn  estuary,  where 
it  does  not  appear  that  there  is  any  hydraulic  control  ac- 
tion at  the  mouth. 

Not  all  vertically  stratified  estuaries  are 
likely  to  be  overmixed.  Rather  obvious  examples  are  Alberni 
Inlet,  which  is  too  deep  to  be  thoroughly  mixed  by  tides, 
and  the  Mississippi  River  passes,  in  which  the  tidal  action 
is  too  small.  However,  in  estuaries  where  there  are  tides 
of  several  feet  or  more,  and  the  depth  is  not  more  then  shout, 
ten  meters,  it  is  reasonable  to  suppose  that  overmixing  may 
occur:  that  i3,  the  salinity  of  the  estuary  upstream  of  a 
transition  is  determined  by  the  control  or  "throttling"  ac- 
tion of  the  transition  on  the  two  layer  flow  through  it,  and 
not  by  the  magnitude  of  large  scale  horizontal  turbulent  ex- 
change in  the  3ense  of  Ketchum  (1951)  and  Arons -Storamel  (1951). 

Hachey  (1939)  has  summarized  the  data  available  on 
the  salinity  of  the  upper  layer  in  St.  Johns  Harbor,  (New 
Brunswick,  Canada).  The  difference  of  salinity  of  the  upper 
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and  lower  layers  is  plotted  as  a function  of  river  discharge 
in  Figure  h jhi.  The  straight  line  is  computed  from  equation 
(4.51.11)  on  the  supposition  that  the  abrupt  widenin'*  at  the 
mouth  of  the  harbor  acts  as  a control,  that  the  width  or  the 
cnannel  is  2,000  feet,  t.ie  depth  of  the  water  is  30  feet,  and 
that  the  mixing  action  of  the  reversing  falls  upstream  is  so 
efficient  that  cvermixing  occurs.  Considering  the  crude  na- 
ture of  both  the  theory  and  data,  the  agreement  is  fairly 
good,  and  the  conclusion  is  that  the  St.  Johns  harbor  is  over- 
mixed. 

Data  on  the  dependence  of  stratification  at  the 

mouth  of  an  estuary  upon  the  river  discharge  is  difficult  to 

find  in  the  literature.  Very  fragmentary  data  suggest  that 

the  following  estuaries  may  be  overmixed: 

St.  Johns  River,  Jacksonville,  Fla. 

Columbia  River,  Oregon 
Savannah  River,  South  Carolina 
New  Waterway,  Holland 
New  York  Harbor,  N.  Y. 

It  is  suggested  that  in  future  surveys  of  estuaries 
where  tne  possibility  of  overmixing  exists,  considerable  de- 
tailed survey  work  be  done  on  sections  likely  to  exert  control 
action.  The  observations  should  give  vertical  salinity  and 
temperature  soundings  over  a complete  tidal  cycle  for  several 
different  stages  of  river  discharge  throughout  the  year.  Con- 
struction of  a graph,  such  as  shown  in  Figure  4.54.1  will  im- 
mediately determine  whether  or  not  the  condition  of  overmixing 
exists  in  a particular  estuary. 
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4 . 6 Salinity  distribution  in  the  Type  2 estuary 

Let  S,  anc  S j.  denote  the  salinity  of  the  top  and 
bottom  layers  respectively;  and  denote  the  discharges 

per  unit  width;  UKr*.  be  the  vertical  entrainment  velocity 
through  the  interface;  and  KL  the  coefficient  of  turbuient 
exchange.  The  equations  of  salt  conservation  may  be  written 

- K Is,  -iv) 

(1) 

5^  + Ks 

The  laws  of  mass  conservation  are  expressed  in  the 
following  relations: 

^ 

ay.  cty  u’ 


It  is  generally  convenient  to  use  ^ as  the  inde- 
pendent variable.,  and  to  write  — K / 

-9 

*'  H< 

dSx 

The  integrals  are  simply 
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where  end  are  constants  of  integration  and  is 

tne  d icons rye  of  the  ri /or  per  er.it  width  of  the  channel. 


fhc  snlir.ity  of  the  ocean  is 


The  discharge  at  the 


moutn  may  be  ‘•V 1'  it  ten  , hence  5 ot  “ 

.touts.  It  is  convenient  to  introduce  ( /^t  and 

or  more  simply  ^ - / - j> , . Thu 


^ 'b 


or  more  simply 


h-sv  . _ /_£n_)7. 
~r~  ‘ l '-p-  / 

*-  (^r-( 


7 (f 


fhus 


In  order  to  illustrate  the  properties  of  these 
solutions,  a numerical  example  is  worked  out.  It  is 
assumed  thut  the  discnarge  of  the  upper  layer  at  the  mouth 
is  ten  times  the  discharge  of  the  fresh  river  water  at  the 
head  of  the  estuary;  the  results  of  the  computation  for 


various  values  of 


are  shown  in  Figure  4.6.1. 


Fig.  4.6.1 
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4. 61  Example : James  River 

An  example  of  the  mixing  processes  in  a Type  2 
estuary  is  obtained  from  rri  tchard  ’ s (1951)  special  study 
of  the  James  River  (see  Section  7. 23).  Figure  4.61.1  shows 
the  mean  salinity-uepth  curves  at  three  stations,  the  es- 
tuary is  clearly  Type  4,  as  shown  in  Figure  7.0.2.  Figure 
4.61.2  shows  the  mean  ebb  and  flood  velocities  for  three 
periods  in  the  James  estuary.  Figure  4.61.3  shows  the 
total  mean  velocity-depth  curves  for  each  period.  If  we  use 
the  customary  right-handed  coordinate  axis,  wi  In  y directed 
downstream,  'j  cross  stream,  upward,  and  denote  the  instan- 
taneous velocity  components  as  “S*  ; the  salt  transfer 

equation  (omitting  molecular  diffusive  transfer)  is: 

It  r 6^-7  a) 

o't  dx  c*}.  Ft 

Each  term  may  be  separated  into  two  terms  of  the 

form:  A = A + A 

where  A is  the  time  mean,  and  A the  instantaneous  devi- 
ation from  the  mean.  Pritchard  (1951)  calls  A a random 
term,  but  this  term  is  more  nearly  a periodic  term.  Tak- 
ing the  time  mean  of  equation  (1)  results  in  the  following 
well  known  form: 


i£  _ - )s  - $£  -u? 
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The  right  hand  member  is  composed  of  three  advection  terms 
and  three  eddy  diffusion  terms..  The  equation  may  also  be 
written  in  the  equivalent  form: 


hr 

dt 


3<r 


(3) 


^In  the  case  of  a Type  1 estuary  the  only  possible  terms 
in  the  right  nand  member  are  the  first  and  fourth,  and  we 
are  reduced  to  tne  case  discussed  in  Sections  4~4  and  4.43). 

Pritchard  has  made  a numerical  study  of  these 
various  terms  in  the  lames  esiuary  to  discover  their  magni- 
tudes, and  takes  only  mean  values  with  respect  to  the  cross- 
stream coordinate  . The  equation  {2}  may  be  written  in 
the  form: 


bir 


■<z  if . jr  2r . jl  (u's  ') 

a*  a* 


2-  fur's  ') 
d± 


(4) 


where  <^1  is  the  breadth  of  the  channel  (a  function  of  £ ). 

The  terms  4i  ,3s /b±  are  computed  from 

a series  of  observations  as  functions  of  ^ and  5L  . The 
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mean  vertical  velocity  i3  computed  from  the  field  of  ZZ 
by  mass  continuity,  and  the  results  shown  in  Figure  4.61.4. 

The  advective  terms  of  equation  (4),  2/Axand  «T  Zf/dt 
may  now  be  computed.  Vertical  integration  of  the  equation  (4) 
permits  evaluation  of  £ (‘ u's  and  then  (ur's')  may  be 


computed. 


Pritchard’s  field  program  was  unuertaken  during  a 

m l^i/dtl  < 10  1 *7ro  . The 


period  when  < 10  °7ro  . The 

term  varied  from  about  5xlO“5  o/oo  sec_l  in 

the  upper  layer,  to  -5xl0“i>  °/ oo^in  the  bottom  layer.  The 
term  Uu.'5‘)/i‘X  was  negligible,  less  than  10“7  °/oo  sec 
The  vertical  advective  term,  ^^~/o  V was  Df  the  order 

of  10“5  °/oo  sec*-^.  

The  vertical  eddy  flux  term  ^ (ur'i  •)  A* 

was 

of  the  order  of  -4xlO"5  °/oo  sec-1  in  the  lower  layer.  The 
vertical  distribution  of  the  vertical  eddy  flux  of  salt 
(ur's')  is  shown  in  Figure  4.61.5.  Since  some  of  the  stations 
were  taken  at  different  phases  of  the  moon,  the  values  of 
eddy  flux  as  functions  of  tidal  velocity  may  be  computed. 
Figure  4.61.6. 

A rough  comparison  with  Section  4.6  may  be  made. 

At  the  halociine  the  value  of  is  about  4xlO~5  ft/sec. 

The  value  of  the  eddy  flux  (ur  's ')  at  the  halociine,  (or 

in  equation  (4.6.1),  is  approximately  40xlO"5  °/oo 
ft/sec.  From  Figure  4.61.1  we  can  see  that  the  value  of 
( 5*1  “ ) is  nearly  -2  °/oo.  Thus,  ^ , as  defined  in  Sec- 
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tion  4.6,  has  a value  of  about  5.  The  salinity  distribution 
for  +j  , in  Figure  4.6.1  resembles  that  of  the  James  es- 

tuary. 

It  is  interesting  to  try  to  compute  the  value  of  UC, 
at  the  halocline  from  Keulegan's  formula  for  mixing.  From 
Figure  4.61.2  we  see  that  u.  - 0.  3 ft.  sec-^-  and  that  this 
is  well  above  Ui-  0.05  ft.  sec  Hence,  by  equation 
(4.8.2)  the  value  of  is  computed  to  be  8x10  ft. 

sec~l,  whereas  the  observed  value,  according  to  Figure  4. 61. 4 
is  more  nearly  4xlO"5  ft.  sec“l. 


**ar  ScM,  V 


Figure  U.61,1  liaan  salinlty-dapth 
curves  »t  thrao  stations 
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Figura  u.6l.<:  iiaan  ebb  voioaitiaa 
(hoavy  linns ) j aeaa  flood  (light) 
for  thro*  periods 


■ - e 


Figure  ii.61.3  it***1  velocitla*  tac 
thrae  periods  la  1*—  Kl»or. 


figure  U.Ol.ii  Abscissa  la 
vortical  velocity  w (with 
nlnus  sipj)  i 1CP  ft  sec-1. 


flCure  U.ul.3  Abscissa 
la  vortical  edly  flux 
TT*  x 10>  o/cc  ft  see-1. 


Fiffav  U.61,6  Uoan  value 
of  vertical  oddy  flux  aa 
function  of  tidal  velocity, 
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4.62  A relation  between  estuaries  with  horizontal  mixing 
and  with  vertical  mixing 


The  salinity  distributions  shown  in  Figures  4.6.1 
and  4.4.1  look  very  similar,  but  come  about  as  the  results 
of  two  entirely  different  models.  The  writer  is  of  the 
opinion  that  ttiere  is  a fundamental  equivalence  between  the 
two  models.  Dr.  Henry  Kiersteed  pointed  out  that  this 
equivalence  may  be  interpreted  in  the  following  way. 

Consider  a two-layer  system  in  which  the  velocity 
of  the  two  layers  are  and  . Suppose  that  the  mean 

velocity  is  /I t . Particles  which  move  randomly  in  the 

vertical  will  have  a downstream  drift  of  . If  the  edciy 
vertical  velocity  of  the  particles  is  U/-  and  they  travel 
on  the  average  the  entire  depth  O^L  , then  a vertical  eddy 
diffusivity  describing  the  flux  of  particles  is  of  the  form 

At  » vrA 

In  moving  up  and  down  the  particles  spend  a time  in  each 
layer  on  tne  average,  so  that  associated  with  vertical  dis- 
placement is  a horizontal  displacement  {u-W’jT  . A hor- 
izontal diffusivity  may  be  expressed  as: 


and  since 


A„  * ('**■>  - L r 

UT  r 2*  0*- 

1 4* 
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The  remarkable  fact  is  that  fa  ^ and  ^^.are  inversely  propor- 
tional. Thus  large  vertical  turbulence  confines  a pollutant 
to  a particular  horizontal  position  which  moves  only  with 
drift  velocity  C . 

Another  way  of  showing  this  odd  equivalence  is  to 
integrate  vertically  the  results  of  Section  4.6,  and  reform- 
ulate them  into  a single  equation  of  the  type  of  Section  4.4. 
For  simplicity's  sake  we  consider  only  the  case,  which  is 
often  realized  in  Type  2 estuaries,  in  which  the  two  layers 
are  of  equal  depth:  D,  * Dj.  » 

The  vertical  average  of  salinity  S'  is  thus 

s = (S,  +Sx)/z 

The  average  velocity  is 


= (g,  +^)/°L 


From  equation  (4.6.4) 


ots 

*ct‘ 


zc  /y*1 

2 C 

L /-/>•  J/'7 


the  quantity 


If  we  write  an  equation  analogous  to  (4.4.8) 


37 


A±  ■ 


t±i± 


i 

! 
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For  large  values  of  4»a  and  O the  equation  is 


simply 


wnich  is  similar  to  the  expression  obtained  by  the  elementary 
particle  analysis  above. 

A r»  important  corollary  of  this  inverse  relation 
between  /\  and  is  that  the  greater  the  tidal  currents, 
the  fresher  the  estuary.  It  should  be  very  interesting  to 
discover  whether  in  fact  an  estuary  of  Type  2 is  freshened 
by  spring  tides,  but  an  estuary  of  Type  1 is  made  more  saline. 
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4.  8 Interfacial  mixing 

The  phenomenon  of  mixing  across  a sharp  interface 
is  peculiar  in  that  it  is  not  a "two  way"  mixing  of  the 
type  described  by  the  usual  theoretical  treatments  of  tur- 
bulence, but  is  a "one  way"  process.  So  far  as  this  writer 
knows,  the  only  study  of  this  type  of  mixing  is  that  by 
Keulegan  (1949),  who  points  out  that  two  forms  of  interfacial 
mixing  may  occur:  "In  one  form,  the  interface  may  be  identi- 

fied as  the  dividing  surface  of  two  layers  of  liquid  with  dif- 
ferent densities,  the  surface  being  one  of  sharp  discontinu- 
ity of  densities  but  not  necessarily  of  velocities.  Ordinar- 
ily the  interface  of  this  type  i3  the  locale  of  internal  waves 
if  the  difference  in  velocities  at  points  on  opposite  sides  of 

the  interface  and  at  some  distance  from  it  is  large.  If  mix- 

ing is  present,  it  is  in  the  form  of  eddies  that  are  period- 
ically ejected  from  the  crests  of  the  waves  into  the  current 

thul  has  the  greater  velocity.  In  the  other  form,  the  inter- 

face is  a layer  of  transition  between  two  currents.  Both  the 
densities  and  the  velocities  change  uniformly  in  the  layer 
that  has  measurable  thickne33.  If  any  mixing  is  present,  it 
is  associated  with  the  momentum  exchange  of  turbulence,  and 
the  regular  pattern  of  internal  waves  is  absent." 

Keulegan 's  study  (1949)  of  mixing  of  the  first  form 
was  made  in  three  flumes  in  which  the  upper  layer  of  fluid  was 
recirculated  many  times  over  the  lower  layer  of  heavier  fluid. 
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I 


The  mechanism  of  mixing  was  the  breaking  of  the  crests  of 
interfacial  waves  into  the  moving  upper  layer. 


i 


I 


\ 

i 

i 

i 


form 


A criterion  of  mixing  has  been  obtained  in  the 


where  subscripts  refer  to  the  two  fluids  in  the  usual  way , 
and  is  the  kinematic  viscosity  of  the  lower  fluid. 

The  value  of  & appears  to  be  about  0.178  in  the 
turbulent  region  ( WT>,/4,>  450). 

Keulegan  defines  the  amount  of  mixing  for  vel- 
ocities above  the  cretical  in  the  same  fashion  as  we  have 
defined  (Section  3.1,  page  6). 

The  law  for  mixing  is  of  the  form 

(U-  > (2) 


"The  bearing  of  the  results  of  the  present  in- 
vestigation on  the  interfacial  mixing  occurring  in  large 
bodies  of  water  in  stratifiod  flow  must  be  discussed. 

"In  the  laboratory  experiments,  where  relatively 


i 


i I 


.->cvji»±o.n  4*  t (rage 


short  flukes  are  used  and  one  of  the  liquids,  the  lower  or 
the  upper,  is  at  rest,  ana  the  other  is  in  motion,  the 
state  of  the  interface  is  one  of  discontinuity  of  density. 

In  the  case  when  the  lower  heavy  liquid  moves,  it  will  be 
supposed  tnst  the  flow  has  continued  for  a long  time  and 
that  the  characteristic  wave  front  is  absent.  Under  these 
conditions,  the  interfacial  stability,  the  critical  vel- 
ocity of  mixing,  and  the  mixing  for  velocities  above  the 
critical,  may  be  studied.  All  these,  however,  refer  to  a 
reach  that  is  to  be  looked  upon  as  a type  of  initial  length. 

"In  the  phenomenon  transpiring  in  natural  environ- 
ments and  thus  involving  large  bodies  of  water  in  stratified 
flow,  it  may  be  assumed  that  conditions  arise  so  that  in- 
itial reaches  are  established.  It  is  expected  that  the 
changes  taking  place  in  the  initial  length  will  be  similar 
to  those  observed  in  a laboratory.  What  is  not  known  def- 
initely in  this  respect  is  the  direct  applicability  oi  the 
results-  to  be  obtained  in  a laboratory  investigation  to  the 
prototype  magnitudes.  Certainly,  however,  a qualitative 
similarity,  at  least,  must  exist. 

"If  that  is  granted,  the  application  of  the  lab- 
oratory results  must  be  restricted  to  a very  short  reach, 
which  will  be  viewed  as  the  initial  reach.  Beyond  this  in 
the  remaining  reach,  which  will  be  of  considerable  length, 
the  conditions  for  mixing  and  the  manner  of  mixing  will  De 
of  a different  type,  obeying  different  laws.  For  the  mix 
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ing  in  the  initial  reach  will  establish  in  the  following 
adjoining  reach  a transition  layer  between  the  liquids,  and 
in  this  transition  layer  the  density  will  vary  continuously. 

As  an  illustration  of  how  this  can  be  Drought  about,  we  may 
visualize  the  following  situation.  A ciirrent  of  fresh  water 
of  depth  H is  flowing  with  a uniform  velocity  over  a pool 
of  heavier  liquid.  In  places  in  the  initial  length,  por- 
tions of  licuid  coming  from  the  lower  pool  will  spread  them- 
selves in  the  upper  current  only  gradually.  Let  it  be  sup- 
posed that  the  spreading  is  proportional  to  time,  this  being 
measured  from  the  instant  of  departure  from  the  lower  licuid. 
The  density  gradient  established  for  this  case  can  be  ob- 
tained readily.  Taking  the  instant  when  the  spreading  is  com- 
pleted and  the  area  covered  is  a square  of  sides  of  length  H , 
the  concentration  along  a vertical  may  be  represented  by 

C r A + 

for  finite  distances  away  from  the  interface. 

"Now,  the  actual  law  of  spreading  for  a given  ac- 
tual case  may  not  be  as  simple  as  in  the  above  illustration. 

As  long  as  it  is  assumed  that  spreading  of  liquids  ejected 
from  below  into  the  upper  current  and  in  the  initial  length 
is  gradual,  a qualitatively  similar  law  for  the  concentra- 
tions as  the  one  mentioned  above  will  be  expected.  But  these 
distributions  emply  the  existence  of  transition  layers.  Where- 
as the  mixing  in  the  case  of  sharp  interfaces  is  brought  about 
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by  tne  ejection  of  eddies  at  the  crests  of  the  internal 
waves,  the  mixing  through  the  transition  layers  must  be 
associated  with  the  momentum  exchange  of  turbulent  motion. 
This  is  a matter  to  be  approached  using  the  basic  ideas  of 
the  Frandtl  and  Richardson  criterion  for  mixing  and  is, 
tnerefore,  a subject-  outside  the  scope  of  the  present  in- 
vestigati on. " 
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4. 9 Measurements  of  microturbulence  in  estuaries 

In  the  theories  presented  in  earlier  sections  of 
this  chapter,  the  turbulent  elements  responsible  for  trans- 
fer of  properties  have  been  considered  to  be  of  a tidal 
period.  If  such  processes  are  said  to  be  associated  with 
"macroturtulence"  then  we  may  call  turbulent  fluctuations 
of  period  much  less  than  a tidal  period  "microturbulence" . 
There  have  been  few  direct  measurements  of  microturbulence 
in  estuaries. 

Thorade  (Rapports  et  Proces-verbaux,  Cons.  Per. 
Intern.  l’Explor  d. 1.  Mer  76,  1931)  discusses  work  of 
Rauschelbach  who  measured  velocity  continuously  in  the  Elbe 
at  about  mid-depth  (total  depth  8m.).  Turbulent  fluctuation 
of  300  second  period,  with  amplitude  of  about  10  cm  sec-1, 
occurred  superposed  on  a mean  velocity  of  about  20  cm  sec  J-. 
Bowden  and  Proudman  (Proc.  Roy.  Soc.  A.  199,  PP.  311-327) 
have  made  an  extensive  set  of  observations  in  the  Mersey, 
and  found  two  periods  of  fluctuation:  (i)  a short  period 
fluctuation  of  a period  of  about  5 seconds;  (ii)  a longer 
period  fluctuation  of  about  90  seconds.  In  both  cases  the 
ratio  of  amplitude  of  the  turbulent  fluctuation  to  the  mean 
velocity  i3  of  the  order  of  0.1  (independent  of  mean  current). 
A positive  picture  of  the  physical  nature  of  the  eddies  (or 
waves?)  responsible  for  these  fluctuations  has  not  evolved, 
nor  do  we  know  tueir  role  in  the  transport  processes.  Much 
observational  work  remains  to  be  done  in  this  direction. 


